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INTRODUCTION 


The  purpose  of  this  study  is  to  develop  an  expression  for  the 


spatial  covariance  for  the  scattered  radiation  when  a  randomly  rough 
surface  acts  as  a  scatterer.  Several  properties  of  the  covariance 
suggest  some  of  the  reasons  for  the  development  of  such  an 


expression.  Some  of  these  properties  include  the  following; 

1.  The  normalized  covariance  is  the  optimum  test  of 

similarity  between  a  pair  of  field  points  by  the  least- 
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squares  criterion. 

2 .  The  covariance  is  the  Fourier  transform  pair  of  the  power 
spectrxim  density. 

3.  It  responds  to  field  fluctuations. 

4.  It  is  less  sensitive  to  inhomogeneities  in  the  propagating 
medium  when  compared  to  the  intensity. 

5.  If  the  covariance  is  only  a  function  of  coordinate 
differences,  then  it  follows  that  there  is  a  uniquely 
determined  fictitious  incoherent  source,  which  if  placed 
at  infinity,  would  produce  the  same  coherence  across  the 
ensonified  scatterer. 

6.  The  correlation  function  obeys  two  Helmholtz  equations 
(two  wave  equations  when  the  disturbance  time  difference 
at  each  receiver  is  larger  than  the  coherence  time)  since 
both  the  field  disturbance  as  well  as  the  correlation 
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between  these  disturbances  propagate  in  the  form  of  waves 
(1.  p.  532). 

It  will  enable  the  computation  of  the  complex  degree  of 
coherence . 


where  the  spatial  covariance  is  -  <p(rj )p*(r2 )> 

and  the  intensity  is  F^^^  ,  i  -  1  or  2,  -  <p(ri)p*(ri)> 

Thus,  this  study  presents  a  useful  tool  in  the  study  of  field 
fluctuations,  source  localization,  and  the  determination  of  some 
physical  properties  of  the  scattering  surface. 

Since  Issakovitch  (2),  Eckart  (3),  and  Rice  (k)  published  their 
formulations  of  the  intensity  of  the  scattered  field  some  35  years 
ago,  a  considerable  number  of  researchers  have  made  improvements  or 
discussed  particulars  of  the  two  basic  methods  that  the  above- 
mentioned  researchers  developed.  The  two  methods  are  the  Tangent - 
Plane  method  (or  Helmholtz-Kirchhoff  diffraction  integral)  and  the 
Method  of  Small  Perturbations.  These  have  surface -roughness  ranges 
of  validity  that  hardly  overlap. 

Several  alternative  models  have  been  proposed,  among  them  the 
Composite-Roughness  Model  (5,6,7)  which  in  essence  merges  the  above- 
mentioned  methods  by  applying  the  Kirchhoff  method  and  the 
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perturbation  method  to  the  range  in  which  they  are  valid.  ;his 
method  assumes  that  the  height  spectrum  can  be  partitioned,  a 
situation  which  turns  out  to  be  possible  in  such  cases  as  the  sea- 
surface  spectrum.  Most  of  the  work  in  this  area  done  by  the 
acoustics  community  concerns  the  sea  surface  or  the  sea  floor  which 
are  assumed  to  have  near-Gaussian  roughness  distributions,  so  that 
most  studies  involve  ideal  pressure -release  or  rigid  Gaussian 
surfaces.  The  problem  of  multiple  scatter,  realistic  boundary 
conditions,  surface  movement,  etc,,  can  make  computations 
exceedingly  difficult  so  it  has  been  customary  to  present  these 
problems  somewhat  isolated  for  the  purpose  of  study. 

At  present,  interest  exists  in  the  development  of  correlation 
and  higher-order  moments  so  that  inhomogeneities  and  fluctuations  of 
the  field  scattered  by  the  rough  surface  can  be  studied.  So  far  in 
the  acoustics  literature,  one  finds  relatively  few  studies  on 
correlation  and  coherence,  most  of  them  being  experimental 
(8,9,10,11,12).  Among  the  theoretical,  Bass  and  Fuks  present  in 
their  monograph  (13)  a  comprehensive  but  general  discussion  on  field 
fluctuations;  Clay  and  Medwin  (14)  develop  approximate  expressions 
for  spatial  and  temporal  correlation  involving  Gaussian  surfaces  and 
include  the  case  of  a  moving  surface  autocorrelation  function. 

Gulin  (15,9,10)  uses  the  Method  of  Small  Perturbations  and  Parkins 
(16)  arrives  at  an  expression  for  spatial  and  temporal  correlation 
using  the  Neumann- Pierson  wave  spectrum.  Kinney  and  Clay  (17)  use 


the  facct-ensemble  method  to  compare  with  laboratory  data  and 
Eckart's  theory. 

Three  methods  are  presented  in  this  study.  The  Method  of  Small 
Perturbations  (MSP) ,  the  Tangent  Plane  Method  (TPM) ,  and  the 
Composite -Roughness  Model  (CRM) .  The  MSP  is  presented  for 
completeness  and  to  emphasize  its  simplicity  and  the  lack  of  the 
need  to  revert  to  the  Kirchhoff  approximation.  MSP  is  useful  when 
the  surface  features  are  small  compared  to  the  source  wavelength. 

It  relies  on  a  small  parameter  expansion  (the  roughness  parameter) 
of  the  surface  to  translate  the  boundary  conditions  on  the  rough 
surface  to  the  (flat)  mean-value  plane. 

The  TPM  is  presented  in  two  different  ways:  the  Standard  Model 
and  the  Slope -Operator  Model.  Both  use  the  Kirchhoff  approximation 
and  have  the  same  validity  range  -  the  small  source  wavelength  limit 
-  and  require  that  the  surface  slope  change  slowly.  The  standard 
model  begins  with  the  Helmholtz  integral  and  relies  on  an 
integration-by-parts  procedure  in  order  to  eliminate  the  dependence 
of  the  resulting  integral  on  the  surface  slope  terms.  The  angular 
dependence  is  expressed  as  a  function  and  is  taken  as  nearly 
constant  so  that  it  can  be  removed  from  the  integral.  This  model 
requires  the  receiver  angles  be  always  very  close  to  the  specular 
direction. 

The  Slope-Operator  model  begins  with  the  first  term  in  the 
iterative  solution  to  the  double- layer  potential  integral.  It  uses 


Gaussian  properties  and  analyticity  in  the  surface  autocorrelation 
to  change  the  order  of  averaging  and  differentiation  in  a  systematic 
way,  so  that  the  end  result  is  an  operator  that  when  applied  to  the 
autocorrelation  provides  the  integral  with  slope  corrections. 

Lastly,  the  CRM  begins  with  the  usual  Helmholtz  integral.  The 
surface  autocorrelation  is  split  into  a  part  corresponding  to  the 
large-scale  surface  contribution  to  which  the  standard  method 
applies,  and  a  part  from  the  small-scale  surface  contribution  on 
which  the  MSP  works  well.  The  problem  is  that  not  all  naturally 
occurring  surfaces  can  be  split  into  two  roughness  regimes.  Another 
problem  had  been  the  choice  of  the  cut-off  surface  wavenumber,  but 
the  implementation  of  diffractive  corrections  circumvented  (7) 
adequately  the  problem  and  improved  the  model's  behavior.  In  this 
study  the  Fraunhofer  phase  approximation  will  be  employed  for  the 
CRM  and  the  diffractive  corrections  will  be  implemented. 

This  study  presents  expressions  for  the  covariance  and 
intensity  using  the  Composite  Roughness  theory,  and  two  formulations 
of  the  Tangent  Plane  Theory.  The  surfaces  in  question  are 
stationary  Gaussian  pressure-release  with  either  Gaussian  or  Pierson 
(18)  ocean  spectrum.  Expressions  for  the  mirror- like  surface 
scatter  are  developed  as  well  so  that  the  coherent  and  incoherent 
components  may  be  studied  separately .  Table  1  provides  the  reader 
with  an  aid  in  understanding  the  organization  of  this  study.  The 
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Table  1.  Diagram  of  the  organization  of  this  study. 

In  parentheses  Is  the  section  In  which  the  covariance 
for  each  model  Is  first  discussed. 
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table  provides  the  equation  numbers  in  parenthesis  which  correspond 
to  each  model.  The  intensity  provides  information  about  the  surface 
elevation,  while  the  covariance  is  more  sensitive  to  surface 
spectrum  and  source  position. 

The  perturbation  approach  is  presented  in  the  first  chapter. 

The  analysis  illustrates  the  fact  that  this  method  applies  to 
surfaces  with  small  crest-to-trough  elevation  amplitudes  but  poses 
no  restrictions  to  the  steepness  of  the  surface  slope.  Further,  it 
supposes  that  the  rough  surface  introduces  "fuzziness”  to  the 
scattered  sound  beam.  This  is  in  contrast  to  the  tangent-plane 
method,  which  proposes  that  the  incident  wave  suffers  Dali-like 
distortion  once  it  impinges  on  the  scatterer. 

A  survey  of  solutions  to  the  Exterior  Oirichlet  Problem  is 
included  in  order  to  emphasize  that  in  TPM  there  are  several 
alternatives  in  the  formulation  of  the  departing  scattering 
integral.  The  Helmholtz -Kirchhoff  integral  lends  itself  to  the 
integration-by-parts  procedure  on  which  the  Standard  Method  is 
based.  The  Double-Layer  Potential  Integral  lends  itself  to  the 
slope -operator  procedure.  The  latter  integral  leads  naturally  to  a 
recipe  for  the  study  of  multiply- scattered  radiation  problems. 

The  survey  suggests  the  TPM  is  a  solution  in  the  geometrical 
acoustics  frequency/roughness  regime.  Further,  the  Kirchhoff 
approximation  is  required,  hence  the  surface  may  have  large  crest- 
to-trough  amplitudes  but  gently  changing  slope.  Lastly,  the  surface 
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must  be  at  least  twice  differentiable  and  the  solution  to  the 
scattering  integral  must  satisfy  Sommerfeld's  radiation  conditions. 

Chapters  4  and  5  illustrate  the  scattering  strength  and 
coherence  predictions  of  the  Tangent- Plane  models  when  a  Gaussian, 
flat  surface,  or  the  ocean  surface  is  involved.  It  is  shown  that 
near  the  specular  direction  the  field  is  mainly  determined  by  the 
coherent  component  and  that  the  incoherent  contribution  becomes 
significant  in  the  field  strength  predictions  away  from  the  specular 
direction.  The  scattering  is  better  predicted  by  the  Slope-Operator 
rather  than  the  Standard  Model  since  the  latter  model  mispredicts 
for  angles  far  away  from  the  specular  direction.  Shown  as  well  is 
the  fact  that  when  the  normal  gradient  d/dn  is  approximated  by  d/dz , 
the  scattering  formula  will  only  display  the  coherent  component  of 
the  field.  This  may  be  adequate  when  the  surface  is  mildly  rough 
and/or  the  autocorrelation  lengths  are  very  large. 

The  Slope-Operator  Model  is  developed  using  the  Fresnel  and 
Fraunhofer  phase  approximations.  Qualitative  comparisons  are  made 
with  regard  to  the  behavior  of  the  model  when  one  or  the  other  phase 
approximation  is  used,  which  show  that  the  Fraunhofer  approximation 
is  adequate  in  the  farfield.  The  plots  for  scattering  strength  show 
that  for  low  grazing  angles  the  Fresnel  expression  consistently 
predicts  lower  energy  amplitudes  as  compared  to  the  Fraunhofer 


expression  and  the  Standard  Model.  The  Fresnel  phase  approximation 
produces  expressions  for  intensity  and  covariance  with  real  and 
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imaginary  parts.  The  author  speculates,  but  never  proves,  chat  the 
imaginary  part  is  less  significant  in  the  farfield  and/or  when  both 
the  autocorrelation  lengths  are  large  and  grazing  angles  small.  The 
coherence  plots  illustrate  that  Che  Standard  Model  and  Che  Slope- 
Operator  Model  predict  similar  results.  The  phase  is  controlled  in 
the  farfield  by  Che  argument  k(r^-r2),  that  is,  the  wavenumber  times 
the  difference  between  the  distances  of  the  receivers  to  ensonified 
surface.  It  is  shown  as  well  that  the  horizontal  coherence  is 
insensitive  to  surface  inhomogeneities  when  the  random  surface 
fluctuations  are  directed  in  the  vertical  direction. 


CHAPTER  I 


PERTURBATION  METHOD  AND  THE  EXTERIOR  DIRICHLET  PROBLEM 


THE  METHOD  OF  SMALL  PERTURBATION  (MSP) 


The  method  of  Small  Perturbation  was  applied  by  Rayleigh  (19) 
to  determine  the  field  scattered  by  periodic  surfaces  and 
subsequently  applied  by  such  authors  as  Rice  (20),  Mandel'shtan  (21) 
to  slightly  rough  random  surfaces.  This  technique,  in  its  simplest 
form,  is  included  here  because  it  forms  a  basis  for  the  Composite- 
Roughness  Model. 

Following  Brekhovskikh  (22)  closely,  the  randomly  rough 
pressure-release  surface  described  by  z-f(r)  with  mean  value 


<r>  -  f  f(r)dS(r)  -  0 
S 


(1.1) 


bounds  the  semi- infinite  region  R.  This  region  contains  a  source 
generating  a  time-harmonic  pressure  field.  For  simplicity,  assume 
the  incident  field  to  be  plane  waves.  Further,  assume  that  the 
total  pressure  at  (r,z)  in  R  satisfies 


Ptotal(r,z,t)-Po(r,z, t)+p(r,z, t) 


z>C  (1.2) 
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where  Pq  is  the  pressure  disturbance  in  the  absence  of  the  rough 
surface.  The  boundary  condition 


Ptotal<*^*r(r)-0 


(1.3) 


is  expanded  in  a  power  series  about  z— 0.  We  retain  only  terms  up  to 
1st  order  (f  is  small): 


P  +  C  ^ 

^total  ^  dz 


total] 


(1.4) 


Substituting  (1.4)  in  (1.2)  and  matching  orders  we  obtain  the 
boundary  condition  at  z-<5‘>: 


0 

1 

0 

(1.5) 

• 

P  -  -  fflCPo) 

(1.6) 

1 

The  solution  to  the  reduced  wave  equation  subject 

to  Dirichlet 

boundary  conditions  is 

(1.7) 

2  2  2 

where  k  -  f  +7 
o  o  o 


^j,»r  -  is  the  horizontal  component  of  kg 

A 

7o*z  -  7(,  is  vertical  component  of  k^ . 
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Substituting  (1.7)  in  (1.6)  and  evaluating  it  z-0 

P(r.0)-2i7of(r)e^^«‘*' .  (1.8) 

Let  A(x)  be  the  surface  spectral  amplitude.  It  is  the  Fourier 
Transform  of  f(r) 

«0 

^(r)  -  S  dxk(x)e^^'^’  .  (1.9) 

Utilizing  this  representation  of  f(r),  we  obtain 

p(r' .z')  -  2i7o  £  (1.10) 

where  7  -  -  ((q+x)^  f  ^  ®nd  Jm{7)>0 

as  the  final  expression  for  the  pressure  at  (r,z). 

2.  Spatial  Covariance  and  Intensity 

The  scattered  pressure  at  (r',z')  is 

p(r',z')  -  2i7c,  J  dx'A(x')e'-^^<’'^^'^‘’'''^'-^'^'  (1.11) 

-CO 

1/2 


with  7'  -  [k2-(^Q+X')2] 

The  covariance  is  the  expectation 

2 

i 

i 

<pp'*>  -  47  //  dxdx'<A(x)A*(x' )> 

■ 

1 

°  ZS' 

« 

\ 

K 

ei^o*(r-r')+i{X- 

(1.12) 
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The  range  of  S  and  2'  are  determined  by  considering  contributions 
from  non- attenuating  waves,  for  which  the  y's  are  real. 

Assume  the  surface  to  be  stationary  so  that  the  surface 
correlation  of  r(r) 

C(rj.)  -  <f(ri)r(r2)>  (1.13) 

depends  only  on  the  separation  rj^  -  r2-ri. 

In  such  a  case  A(x)  and  A*(x' )  are  delta  correlated  (13,  p.  43) 

<A(x)A*(x')>  -  W(x)5(x*x') .  (1.14) 

where 

OO 

W(x)  -  X  dr^  C(ri)e'^^”'^  (1.15) 

is  the  spatial  spectrum  of  roughness. 

Substituting  (1.14)  in  (1.12)  and  integrating  with  respect 

to  x'  . 


<pp'*>  -  47  /  dxW(x)e^^^'‘’^ 

2 


i(X+Co) •ti+i(7z-7'z' ) 


(1.16) 


is  the  covariance.  When  (r,z)— (r' ,z' ) ,  the  intensity 


<pp*>  -  47  f  dxW(x) 


(1.17) 


is  obtained. 


-«r\j  iTU  V  j  ir-j  ir  -  ^ 
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3 .  Discussicn 

Equation  (1.10)  implies  that  to  the  first  approximation  each 
Fourier  component  of  roughness  acts  as  a  virtual  source  giving  rise 
to  its  own  scattered  plane  wave  propagating  in  the  direction  (^,r), 
which  obeys  Bragg's  Law  (22,  p.  177): 

€  -  Co  +  X  •  (1-18) 


Equation  (1.16)  and  (1.17)  show  that  the  scatter  is  mainly 
determined  by  the  surface  spectrum,  not  the  height  distribution,  in 
the  regime  in  which  this  method  applies.  The  main  feature  of  MSP  is 
the  translation  of  the  boundary  condition  on  z-C(r)  to  z-<C>-0  by 
means  of  a  Taylor  series  expansion.  This  fact  restricts  its 
applicability  to  wavelengths  that  are  large  compared  to  the 
correlation  length  p  of  the  surface  (13,  p.  82),  or  more  precisely, 
to  the  range 


R  -  2koasinSg  «  1 


(1.19) 


where  R 


is  the  Rayleigh  parameter 

is  the  source  wavenumber 

is  the  rras  height  of  the  surface 

is  the  grazing  angle  of  the  incident  wave. 


15 


Physically  this  amounts  to  each  irregularity  weakly  distorting 
the  field  incident  on  it,  and  that  in  the  frequency  range  of 
validity  the  roughness  of  the  surface  is  shallow  as  seen  by  the 
incident  wave. 

Two  alternatives  to  the  MSP  will  be  presented  next,  one  being 
the  popular  "Tangent  Plane  Method"  (TPM) ,  which  assumes  that 
reflection  takes  place  at  each  point  on  the  surface  according  to 
geometric  acoustics  (the  incoming  wave  sees  a  surface  chat  very 
gradually  changes  and  reflects  from  Che  tangent  plane  at  each  point 
of  Che  surface) .  The  problem  with  this  method  is  that  the  Kirchhoff 
approximation  is  needed  to  predict  the  field  at  the  boundary.  The 
other  model  is  a  combination  of  the  MSP  and  Che  Tangent  Plane  Method 
and  it  assumes  that  Che  surface  is  to  be  considered  as  the 
superposition  of  gently  changing  large-scale  roughness  with  an 
overriding  small-scale  ripple. 

The  MSP  has  been  presented,  for  simplicity,  with  an  incident 
plane  wave.  However,  the  source  could  just  as  well  generate 
spherical  waves  and  have  a  beam  pattern.  In  both  the  Tangent  Plane 
method  and  the  Composite-Roughness  Models,  source  characterization 
will  play  an  important  role  in  the  following  presentation,  but  this 
need  not  be  so  in  the  general  case  (23) . 

Phenomenologically,  it  is  interesting  to  note  that  the  MSP  and 
the  TPM  differ  in  their  prediction  of  the  effect  the  scattering 


surface  has  on  the  incoming  wave:  MSP  predicts  that  the  inci'^ent 
wave  becomes  "fuzzy"  or  diffuse  on  reflection  whereas  TPM  predicts  a 
scattered  field  as  a  superposition  of  sharp  images  of  the  incident 
wave  producing  a  "distorted"  but  well  outlined  incident  wave.  The 
composite-roughness  result  is  a  combination  of  both  effects. 

B.  THE  EXTERIOR  DIRICHLET  PROBLEM 

1 .  Statement  of  the  Problem 

The  total  acoustic  pressure  at  r  in  a  semi- infinite  region  K  of 
isotropic  media  will  be  assumed  to  be  the  point-wise  superposition 

Ptotal(’^'t)-Po(r,t)+p(r,t)  ,  (1.20) 

where  Pq  and  p  are  the  contributions  from  the  source  at  r'  and  the 
boundary  respectively  (see  Fig.  1).  Our  attention  ce .iters  on  the 
scattered  field  p,  which  in  the  time-harmonic  case  obeys 

(A+k2)p(r)e'i-“*^  -  0  in  R  (1.21) 

the  reduced  wave  equation,  with  Dirichlet  boundary  condition  on  the 
randomly  rough  interface 


p(r")-0 


r"  G  S 


(1.22) 
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and  Sommerfeld' s  Radiation  conditions  (24,  p,  429) 


The  double- layer  potential  with  continuous  density  ^  is  a 
solution  to  the  Exterior  Dirichlet  Problem  (1.21)  (1.22),  provided  0 
is  a  solution  of  the  integral  equation 

A 

2p(r)  =  0(r)+2  J  i/)(r" )n"  .V"0(r ,  r")dS"  when  r  e  X  (1.25) 

2  and  r"  ^  r  . 

This  equation,  in  turn,  has  a  solution  if  the  surface  belongs 
to  a  class  of  Lyapunov  surfaces  (27)  and  if  the  density  0  is  bounded 
everywhere  on  2.  One  important  implication  of  this  restriction  is 
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> 

1 

( 

that  the  correlation  function  must  be  analytic,  and  secondly  the 
function  f (r- ■ )  describing  the  surface  must  be  Holder  continuous 
(28). 

Both  p  and  ^  are  unknown  on  the  surface ,  but  an  approximate 
expression  for  the  pressure  for  points  on  the  surface  is 

p(r)  -  -Po(r)  r  e  X.  (1.26) 

which  is  commonly  known  as  the  geometrical  optics  approximation  (or 
one  of  the  Kirchhoff  approximations)  because  it  is  only 
approximately  true  for  k|r-r"|»l  --  the  high  frequency  regime 
The  physical  intepretation  of  (1.26)  is  chat  at  points  corresponding 
to  the  surface,  the  field  equals  the  incident  pressure  at  those 
points  in  the  absence  of  the  surface.  Further,  when  r  is  on  the 
surface  S” ,  but  not  in  2,  the  pressure  is  zero  at  r.  Thus 
specifying  that  shadows  be  infinitely  sharp. 

Substituting  (1.26)  in  (1.25)  and  performing  the  iteration 

0(rG2)  -  V’o  +  ^1  +  ^>2  .  .  .  V’m  ■  •  ■  . 

where  V’o(J^)  ”  '  2  Po(reZ) 

V>i(r)  -  -  2  /  ^o(ri)ni*Vi,^(r,ri)dSi 

Si 

V>n,(r)  -  -  2  f  V*  ,  (r  )n  ‘V  (^(r,r  )d 
'  i,  ^m-1  m  m  m  m 

2^ 


res 


r,e  Si 


IS 


re  S„ 


(1.27) 


Retrining  only  the  first  two  terms  in  (1.27)  and  substituting  in 
(1.24) 


p(r)  -  *  2  /  Po(r")n"*V’.^(r,r'')dS" 

4  /  {  /  Po(r")ni.v"^(r,ri)dSi”}  n" .7"^(r, r" )dS"  -  ...  (1.28) 

2  '■  S  ^ 


The  first  order  contribution  is  recognized  as  the  Rayleigh 
Diffraction  Integral.  Equation  (1.28)  states  that  the  field  at  r  is 
a  linear  combination  of  n-order  multiple  scatter  terms.  The  first 
term  represents  the  singly  scattered  contribution  since  the 
integration  over  2  means  that  the  contribution  of  every  point  r"  is 
summed  coherently  to  yield  the  first-order  contribution.  The 
second,  that  part  of  the  disturbance  at  points  r"  in  2^ ,  reradiates 
in  the  direction  of  r"  in  2,  so  that  the  second  order  contribution 
to  the  pressure  at  r  is  the  integral  over  2  of  the  coherent  sum  in  2 
of  the  reradiation  to  every  point  r" .  Higher  order  terms  in  (1.28) 
represent  even  more  complicated  reradiation  contributions.  This 
method  is  particularly  well  suited  for  the  study  of  multiple 
scattering,  especially  since  iteration  and  superposition  are  in 
principle  relatively  easy  to  implement  on  a  computer. 

In  principle  it  is  also  possible  to  formulate  the  solution  to 
(1.21  and  1.22)  by  the  use  of  a  Single  Layer  (29)  or  Simple  Source 
Distribution  Integral 
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p(r)  -  /  0(r"><^s(r.r‘')dS"  in  »,  r"  e  S".  (1.29) 

Z 

Equation  (1.29)  is  a  solution  to  the  Exterior  Dirichlet  Problem 
in  question,  provided  ^  is  a  solution  to  the  homogeneous  integral 
equation 


/  0(r")^s(r.r")dS"  -  0  .  (1.30) 

Z 

The  task  would  be  to  find  4s  which  is  identically  zero  at  the 
randomly  rough  boundary  and  that  satisfies  radiation  conditions. 

This  makes  this  formulation  very  unattractive  and  will  not  be 
discussed  further. 

Finally,  we  address  the  Helmholtz  Integral  formulation  for  the 
scattered  field 

p(r)  -  /  {  0(r,r")n"«V'’p(r")  -p(r")n"*V"(^(r,r")|  dS"  r"  e  Z 

(1.31) 

This  is  the  most  frequent  starting  point  for  authors  discussing 
the  Tangent  Plane  Method.  The  formula  may  be  derived  using  the 
divergence  theorem  and  Green's  identities  (30,  p.  803;  24,  p.  40). 

Equation  (1.31)  requires  that  the  pressure  and  its  normal 
derivative  at  each  point  of  the  surface  be  known.  For  a  randomly 
rough  surface  this  can  only  be  determined  approximately,  for 
example,  invoking  the  Kirchhoff  approximations. 
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P(r")  -  -  Po(r")  (1.32) 

^  A 

n".V”p(r")  -  2n"*7"Po(r")  (1.33) 

part  of  which  was  used  in  connection  with  the  dipole  layer 
formulation.  If  these  approximations  are  made,  (1.31)  is 
transformed  into 

-  2  ;  (Po^)dS".  (1.34) 

Pq  and  4  are  the  doubly  differentiable  incident  pressure  and  the 
Fundamental  Green's  Function  respectively,  n"  is  the  outward  normal 
of  the  randomly  rough  surface.  This  formulation  could  incorporate 
multiple  scattering  if  the  full  integral  equation  that  leads  to  the 
Kirchhoff  approximations  were  implemented  (29) . 

Alternatively,  can  be  chosen  to  be  the  solution  of 

(V2  +  k2)^^(r,r")  -  -  5(r-r")  (1.35) 

which  vanishes  at  the  flat  surface  z-0.  Under  this  condition  (1.31) 
is 

P(r)  -  -  X  P(r")f-  0vi(r.r")  IdS"  (1.36) 

S  z-0 

The  solution  of  (1.36)  is  arrived  at  by  first  making  a  Taylor 
series  expansion  for  small  a,  the  rms  height  of  the  surface,  much 
like  in  MSP,  that  would  translate  the  boundary  condition  on  the 
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rough  surface  back  to  the  surface  z-0.  This  technique  is  limited  to 
wavelengths  much  greater  than  a  but  does  not  rely  on  the  Kirchhoff 
approximations . 


The  previous  section  entunerates  some  of  the  less  esoteric 
formulations  of  the  solution  to  the  Exterior  Dirichlet  Problem. 

Each  one  of  these  is  used  as  a  starting  point  for  the  development  of 
different  formulations  of  the  TPM.  Equation  (1.34)  and  the  first 
term  of  (1.28)  (it  is  assumed  that  multiple  scatter  does  not 
significantly  perturb  the  scattered  field)  will  be  developed  further 
in  the  next  chapter.  They  both  incorporate  the  Kirchhoff 
approximations,  but  differ  in  their  starting  points  as  well  as  in 

A 

the  evaluation  of  the  quantity  n  •  k  and  the  phase  in  the 
exponential.  Since  Kirchhoff’ s  approximations  are  invoked, 
diffraction  shall  be  ignored.  Thus,  surfaces  to  which  the  Tangent 
Plane  Method  applies  may  have  large  crest  to  trough  differences ,  but 
slope  changes  must  occur  very  gradually.  Furthermore,  it  is  assumed 
that  the  source  sees  every  point  on  S.  Otherwise,  eqs .  (1.32)  and 
(1.26)  are  replaced  by 


p(r")  -  -S(r")Po(r") 


(1.37) 


where  S  is  a  shadowing  conditional  probability  function  (31)  that 
depends  on  the  source  depression  angle  as  well  as  the  height  and 


slopes  of  the  surface.  It  is  clear  that  this  factor  becomes  more 
important  for  smaller  source  depression  angles,  but  in  this  work  we 
shall  assume  it  is  always  equal  to  unity. 


CHAPTER  II 


THE  TANGENT  PLANE  METHOD 


A.  THE  TANGENT  PLANE  METHOD 


The  starting  point  is  the  Rayleigh  diffraction  integral  (1.28) 
and  the  Helmholtz  Integral  result  (1.34) 


p(ri)  -  -  2  /  Po(r")n".7-,^(ri.r")dS" 
2 


P(ri)  -  2  J  ^Po(r")4S(ri,r")j  dS"  . 


(2.1) 


(2.2) 


As  illustrated  in  Figure  1,  the  source  at  r'  generates 
spherically  spreading  waves,  limited  by  the  beam  function  D" : 


(2.3) 


where  R'  -  |r'+r"| 


The  Fundamental  Green's  Function  is 


0(ri ,r")  - 


4^|ri-r"| 


^ik|ri-r"|  _  1  ikRi 


(2.4) 


and  has  a  gradient  with  normal  projection 


n"  •  7"(^(r,r")  -  in  • 


(^-  ifej  - 


(2.5) 


Assume  that  kR^  »  1  so  that  expression  (2.5)  can  be  approximated 


Further,  the  source  is  assumed  to  have  a  narrow  beamwidth  so  that 
the  approximations 

R'  -  r' 

-  r^  are  valid.  (2.9) 

Another  convenient  change  that  can  be  made  on  (2.1)  by  this 
choice  of  D  is  to  change  the  limits  of  integration  so  that  the 
integration  is  perfoirmed  over  the  whole  surface.  It  is  this 
assumption  that  makes  it  evident  that  the  choice  of  source 
beampattem  plays  an  important  role  in  the  development  presented 
here. 

3.  Fraunhofer  and  Fresnel  Phase  Approximations 

Both  and  Pq  (ignoring  the  source  function  D)  in  (2.3)  and 
(2.4)  represent  spherically  spreading  waves.  The  phase  in  (2.3)  and 
(2.4)  is  found  to  be  (see  Fig.  2) 

ikR'  -  ik[x"2  +  (r'cosS+y")2  +  (r ' sinS -f" )2 ] ^ ^ 

ikR  -  ik[ri^cos5 j^cos^^ -x")^  +  (ticos^j^sini^^ -y")^  + 

(riSini?i-?-")2]^'^  (2.10) 

These  can  be  approximated,  retaining  the  first  two  terms  of  the 


binomial  expansion,  by 


IL'"  XT  X,"  Ttr  XTT  X.*  icr  wj*  KT.  ILI  XT  V-  RJT  •J'  KT. 


ikR'  -  ik|r'+y''costf-{‘"sintf  '*’'^7  (x"^+y"^)j 


(2.11) 


ikR  -  ik^ti  -  (costfiCos^jX"+costfiSin^iy+sin5if" )  +  ^^(x"2+y"  2)j  . 


(2.12) 


In  the  future  these  two  phase  appear  added 


ik(R'+R)  -  ik[ri+r‘+aix"+)3iy’’-7if"]  +  ikPj  (x'’2+y''2 ) ,  (2.13) 


where  -  cos^^cos^^ 


5i  -  -  costfjsin^i 


4i,  -  (cos5-5j) 


7i  -  sintf  +  sin^i 
r,  +r' 


This  expansion  is  valid  when  r'*^  (x"^+y"^+2r' (y"cosS-{'"sinS) } 

<  1  and  r^^  {x"*+y"2-2r^  (e^x"+5^y"+sin^5'") )  <  1,  which  physically 
implies  that  the  source  and  receivers  are  far  enough  so  that  the 
largest  dimension  of  the  insonified  area  is  still  smaller  than  the 
source/receiver- to-surface  distance.  An  alternative  to  the  binomial 
expansion  is  the  replacement  of  R+R'  by  a  power  series  in  r"  about 
the  origin  (33,34).  The  binomial  expansion  is  admittedly  poorer 
than  the  series  expansion,  but  it  is  considerably  easier  to 
implement . 

When  the  source  is  far  away  so  that  kR'  »  1  then  the  wave 
fronts  of  Pq  are  locally  nearly  plane  and  Pq  can  be  approximated 


■r.  •r.  .V  -r. 
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Po  - 


(2.14) 


r' 


with  k'  -  k  ”, 
r ' 


Similarly  ^  can  be  approximated  if  kR]^  »  1  so  that 


J - iki-(r,-r") 


4xR, 


(2.15) 


When  the  phases  of  (2.14)  and  (2.15)  are  combined,  the  term  in 
square  brackets  in  (2.13)  is  obtained.  This  is  labeled  the 
Fraunhofer  phase  approximation  and  physically  suggests  the 
assumptions  made  on  the  shape  of  the  wavefronts.  When  the 
approximations  in  (2.14)  and  (2.15)  are  not  made,  i.e.,  the 
wavefronts  are  allowed  to  have  some  curvature,  then  at  least  second 
order  terms  must  be  included.  When  second  order  terms  are  retained, 
as  in  (2.13),  the  Fresnel  phase  approximation  is  being  made.  This 
approximation  somewhat  complicates  computations  of  the  moments  of 
the  pressure  and  shall  only  be  given  some  attention  here  when  the 
Slope -Operator  Model  is  discussed. 

Before  proceeding,  (2.1)  and  (2.2)  are  restated  incorporating 
the  changes  discussed  so  far: 


ik(r, +r' ) 
e _ ^ 


P(ri)  -  27ir  r'  dS"n"-kiDe 

^  -  eo 


ik(aj^x"+i9iy"  -7i  f"+Pi  [x"^+y"^  ] ) 


•  00 


(2.16) 


31 


.ik(ri+r') 


P(r,)  -  (2.17) 

1  .(B 


with  n" 


-3r".  "  ^  " 

(3x"^  ^  <3v"^  y-^ 


at"  2 


ar“  * 
(^^)  + 
^dy"^ 


1] 


1/2 


and  dS"  -  |Vf"|dx"dy"  . 

Equation  (2.16)  corresponds  to  the  Rayleigh  Diffraction  Integral 
expressed  so  that  either  Fraunhofer  or  Fresnel  phase  approximations 
can  be  considered.  Equation  (2.17)  is  the  Helmholtz  integral  result 
and  has  the  Fraunhofer  phase  approximation. 

The  evaluation  of  (2.1  and  2.2)  shall  be  performed  in  two 
different  ways.  One  approach  shall  be  labeled  "Standard  Model"  and 
starts  with  (2.17).  The  "Slope  Operator  Model"  uses  (2.16)  first, 
with  Fraunhofer,  and  later  with  the  Fresnel  phase  approximation. 

The  latter  model  does  not  make  a  large  grazing  angle  approximation 
as  will  be  the  need  in  the  Standard  Model  in  order  to  permit  the 

A 

evaluation  involving  the  normal  n" . 


B.  DERIVATION  OF  THE  INSTANTANEOUS  AND  MEAN  SCATTERED  PRESSURE 
1 .  Standard  Method 
a.  Development 

This  approach  has  been  presented  with  slight  variants  a  great 
many  times.  Among  the  earliest  researchers  to  apply  it  to  the 
computation  of  intensity  are  Issakovitch  (2) ,  Eckart  (3) ,  and 


>/■> 


-’vS 


'53 


r\'S 

li  •  ■  " 


•  ■  ■« 


>>1 


Brekhovskikli  (35) .  Here  we  follow  the  approach  presented  by  Tolstoy 
and  Clay  (32)  and  use  the  geometry  of  Figure  2. 


Using  the  Fraunhofer  approximation,  the  difference  in  the  wave- 
vectors  is 


2^  ^  A  A 

^  (k'-kj)  -  -  £jX  +  (cosfl-fij)y  -  -y^z 

with  -  costfj^cos^^ 

Sj  -  costf^sin<^j 
7i  -  sine  +  sin^i 

so  that 

ik<r'-vrj^>  <»  «. 

P(ri)  -  -  k  /  dy"  /  dx"D" 

^  .A  *  A 


(2.18) 


(«ir"x+^ir 


(2.19) 


with  -  -  £j 

Pi  -  (cose-5j^) 

»t  II 

Cjj,  Cy  are  respectively  partials  with  respect  to  x 
and  y"  of  C"  . 

The  first  two  terms  can  be  simplified  using  integration  by  parts. 
The  first  term  gives 


k  /  dx"a,D"e"^i^"(e-"‘'^i«‘f”dx«)  -  -  .-  oe"  ^  ' 

L  '  ^  ^T-i 


+  k  f  dx-D-  -  +  f  dx"  ?-  3ik(a,x"-7,r") 

.00  '>'1  .®  ^'>'1 

(2.20) 


Substituting  the  explicit  form  of  and  D"  gives  the  following 
result: 


k  /  dx"aiD"e'^^i^"(e'^^”’'i^"c”dx)  -  0 


<*>  D"a  .,  .  ®  dx"x"D"a  .,  ,  ,, 

+  k  ;  dx - i-  e  ^k(a,x«-7if').2k  /  - r-TT"^  ^ik(a,x" -7: f " ) 


-  F(U)+G(U) 


For  large  kU,  G(U)  falls  off  faster  than  F(U) .  By  Schwarz 


(2.21) 


inequality 


IG(U) I  _1_ 
|F(U)|  “  kU 


so  that  G(U)  is  negligible  if  kU  »  1  is  enforced.  The  procedure 
is  used  to  perform  the  y"  integration,  yielding 

ik(r'+ri)  ca  ®  (a"  +  +  7?) 

*  ~27r'7i —  ^  ^  - 7i - ^ 

-  OD  -00  ^  ^ 


(2.22) 


provided  that  kV  »  1  as  well. 
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Finally,  substituting  (2.18)  in  (2.22)  and  imposing  the 
restriction  that  changes  in  the  incident  and  receiver  grazing  angles 
be  restricted  to  a  small  range  so  that  the  angular  dependence  may  be 
removed  from  the  integral,  the  final  result  is 


ik(r'+ri  )  ® 

~  2^  r>  '  ^dx"dy"  , 

1  -ca 

(2.23) 


(2.23)  is  valid  for  angles  in  the  neighborhood  of  the  specular 
direction  and  kU,kV  »  1. 

If  the  surface  is  flat  f"-0  and  the  scattered  pressure  from  a 
flat  pressure -release  surface  or  a  "mirror"  is  obtained: 


"  '^^2x!  r'  '  D"e^^^“i'""''^iy"^dx"dy"  (2.25) 


b .  The  Mean  Scattered  Pressure 

The  surfaces  to  be  dealt  with  will  be  random  stationary, 
described  by  C"(x",y").  C’"(x",y")  is  a  sample  of  a  stochastic  process 
and  the  scattered  pressure  is  therefore  a  random  quantity.  The 
ensemble  average  or  "mean"  pressure  is 


(2.26) 


The  interchange  between  the  expectation  and  integrations  is  possible 
when  there  is  uniform  convergence.  The  expression  <e  >  is 

recognized  as  the  characteristic  function  of  the  random  quantity 
r"(x".y") 


(2.27) 


where  a)(f")  is  the  probability  density  function  of  heights  f".  As  a 
special  case,  consider  the  Gaussian  probability  function 


w(C") 


-C"*/2o2 


Jlira^ 


(2.28) 


Substituting  (2.28)  in  (2.27)  we  obtain 

2 

where  g  is  the  roughness  parameter,  defined(*)  as  g  - 


(2.29) 


The  mean  scattered  pressure  for  a  Gaussian  distributed  surface  is 


<p(ri)>  - 


(2.30) 


Raadars  familiar  with  acoustica  litaratura  that  deals  with 

2  2 

scattering  Croto  rough  surfaces  «#1I1  recognize  the  quantity  k  7  ^ 

2  ^ 
as  defined  by  g  .  In  this  study  the  square  is  omitted  so  that  g 

*  k^7^(7^  . 

1 


thus 
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.  g/2 

The  factor  e  ^  is  proportional  to  the  degree  of  coherence  in  the 
scattered  radiation.  It  tends  to  one  for  mirror  surfaces  and  to  zero 
as  the  roughness  increases.  It  is  worth  noting  that  (2.30)  implies 
that  the  average  pressure  is  independent  of  the  spectral  composition 
of  the  surface. 


The  slope  operator  itself  will  be  presented  in  the  context  of 
higher  order  computations.  Described  here  is  the  derivation  of  the 
field  using  Rayleigh's  diffraction  formula.  The  result  presented  here 
will  be  used  later  when  the  operator  is  derived.  This  approach 
appears  a  few  times  in  the  acoustics  literature,  notably  in  Welton's 
report  (25),  as  well  as  in  Boyd  and  Deavenport’s  paper  (36). 
a.  Development 

The  dot  product  in  eq.  (2.16)  has  components 


(n"»ki)i  -  —  (riti-x")?- 


(2.31) 


(n".ki)j  -  ^(ri5i-y")r^ 


(2.32) 


(n"*ki)z  -  — (r^sin^i -C”)(-l) 


(2.33) 


Let  ■  aj^x''+/9iy'' -7^  (■"+?!  [x"*+y"*  ] 
Equation  (2.16)  is  therefore 
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ik(ri+r'),  «>  ® 


2^ir,r'  J'  ^ 

i  .00  -  00  A 


X  {(ri«i  -x")f  +(ri5i  -y'’)r„-(riSintfi  -f") ) 
X  y 


(2.34) 


Rearranging 


P(ri)  - 


rir' 

k 


//  dx"dy” 


(«ir  +i5if„-sin«i) 


X  ^'y 


It  It 

j-r  M 


rT  .  „u  „  D1  ik*"  r 


(2.35) 


The  second  term  can  be  ignored  if  the  beam  is  sufficiently  narrow  so 


that  x"«ri  and  y"«ri  .  Further,  note  that  this  term  is  of  order  r^' 


smaller  than  the  first  term.  Thus, 


^ik(ri+r'),  <» 


(ri)  -  j.< —  //  D"e  (eif^+5  -sin^i)dx"dy’’.  (2.36) 

1  .00  ^ 


The  pressure  reflected  by  a  plane  pressure-release  surface  is 


obtained  by  letting  f  approach  zero: 


_  ike^^^*^^'*'^'^in^.  p7  ik(aiX"+^iy"+Pi[x"2+y 

2xr,r'  JJ  ® 

^  .m 


^  1 ) 

‘  '^dx"dy" 


(2.37) 


In  the  farfield  the  Fraunhofer  approximation  applies,  thus 


[x"^+y"2 ]-0.  Attention  to  the  matter  of  phase  approximations  shall 


be  left  to  later  sections. 


ElWCSWiTOWt 


Again,  assume  uniform  convergence  so  that  the  expectation  and 


I 

spatial  integration  operations  may  be  exchanged.  Thus, 


1  .«Q 


X 


<  e 


-ik7ir 


li£x 
sintf , 


n 


sinS  j 


) 


>dx"dy" 


(2.38) 


When  C"  is  statistically  stationary  then  the  average  can  be 
separated  from  the  integrations  in  x"  and  y"  (25) 


<P(r.,>  -  )  >P„  <2.3„ 

For  the  special  case  of  Gaussian  random  process  with  zero  mean 
height  and  slopes,  (2.39)  yields 

<p(ri)>  -  (2.40) 

as  before. 

3.  Discussion 

A 

So  far  the  matter  of  evaluating  the  dot  product  involving  n" 
has  been  illustrated.  The  results  in  sections  1  and  2  differ  in  the 
evaluation  of  the  dot  product  because  the  integrals  were  different  to 
begin  with.  These  discrepancies  are  apparent  when  comparing  (2.26) 
and  (2.38)  as  well  as  the  plane  pressure  release  results. 


Some  authors  choose  to  approximate  the  average  of  the  random 

A  A 

vector  n"  by  the  z  coordinate  unit  vector  (3)  when  the  slopes  are 
gently  changing  and  have  zero  mean  value  on  surfaces  for  which  the 
Tangent  Plane  Method  applies.  The  slope  corrections  are  increasingly 
important  at  lower  grazing  angles,  far  from  the  specular  direction, 
and/or  when  the  height  surface  spectrum  content  is  high.  In  the 
Standard  approach  the  approximation  is  reasonable  since  the  model  only 
applies  to  the  neighborhood  of  the  specular  direction,  but  it  is  a 
severe  restriction  for  the  Slope-Operator  Model  since  it  is  valid  for 
a  larger  range  of  grazing  angles  and  incorporates  slope  corrections. 


in 


■Vl 


^  •  v' 


ft 


iHAU 
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CHAPTER  III 


SPATIAL  COVARIANCE  AND  INTENSITY 


The  general  expressions  for  covariance  <p(r]^)p*(r2)>  and 
intensity  <p(r)p*(r)>  for  the  scattered  pressure  involving  surfaces 
with  stationary  Gaussian  height  distributions  will  be  presented  using 
the  Standard  Slope-Operator  and  Composite-Roughness  models. 

A.  STANDARD  MODEL 
1 .  Formulation 

Using  (2.23)  with  (2.24),  the  pressure  at  r2  is 


,,  ik(r'+r2 )  «  -i  / 

2  .«o 

(3.1) 


fjCe)  - 


l+sintfsinSj  -costfcos^2  sin(^2 
sinS+sinSj 


(3.2) 


Hence ,  the  spatial  covariance  is 

k2fj(e)f2(e)e^^^’^i'’^2^  « 

<P(ri)P*(r2)>-<  ■  4x2r'2rjr2 - 


^ik(aiX'  ■  -a2X'  •  ■+0iy  •  -fi2y  ’  ’  ).g-ik(7ir'  '  '  '  ) 


dx'  ’  dy  ’  dx'  ’  ■  dy  ’  ■  > 


As  usual  <  >  means  ensemble  average  and  *  means  conjugate.  Fc 
brevity  the  symbol  rj^2  shall  replace  <p(ri  )p*(r2  )> ,  the  covariance, 
and  or  r2  2  for  the  intensity  at  and  r2  respectively. 
Expediency  in  computations  is  improved  by  changing  coordinates  to 
center-of-mass  and  difference  coordinates.  Let 


and 


Thus 


y 


V  ■  +V' ' 
2 


X  -  X' ' -X’ ' '  y  ”  y  '  *y’ ' ' 

A  A 

rjL  -  xi+yj 

A  A 

Rx  -  xi+yj 


Since 


a(x'  '  ,X'  ‘  •  ) 

a(x,x) 


1  and 


a(y, .  ^  y  •  ■  ) 

a(y.y) 


1  then 


dx'  ■  dy  ■  dx'  ■  •  dy  •  ■  -  dxdxdydy 


Let 


dr^  ■  dxdy 
dRx  ■  dxdy 

9 


9 


Finally,  let  a  - 


SO  that  s  ■  3L+pj 


s  ■  ofi+^j 


Assuming  uniform  convergence,  the  expectation  is  taken  inside  the 
integral  in  (3.3).  We  adopt  the  changes  in  coordinates  and  notation 
to  get 


(e)f2  ®  CO  ^  _ 

Tu  -  - ^  <  X  c  r,  s) 


.00  -CO 


■ik(7if’ ■ -yz^'  ‘ '  )> 


(3.4) 


where 


D  -  exp{-2(^  + 

U2  v* 


1  v2 

D  -  exp{  -  j[ —  +  ) 

1)2  V2 


The  expression  in  angular  brackets  is  recognized  as  the 
characteristic  function  of  the  bivariate  probability  distribution 
function  w  of  the  surface,  defined  as 


^-ik(7ir”-72r-')>  _ 


-00  -CO 


(3.5) 


Since  it  is  assumed  that  the  surface  is  (.lussian  in  distribution, 
homogeneous  and  stationary,  the  characteristic  is  computed  by  first 
expanding  (3.5)  in  a  Maclaurin  series: 


<e‘i^(^i^'"'^2^’")>-<l-i(  )-  ^[7^"=“-27i72f"r--+72r'”"]+  i^( 


3! 


IT  t'l'iy'  ’  ■^T'ly  ’  3^2?'  ’  ’+67iy'  •  ^72r'  ■  ■  ^-^7i!:'  ■  727'  ’ '  +737'  ’ 


(3.6) 


Following  this,  the  expectation  is  computed  term  by  term: 

T'af  )v,  _  ^  27^72^^0+720^) 


+  ^  (37  (T*'-127  7^0  C+I87  7^0^0^-127,7  0^0+37  o*  )  .  .  . 
4!  '1  i.  12  4  '1  2  4  '^  2  4  14 


—  l-k^o^ 


r  2  2  \  r  Z  Z  'v 

K‘^^2  k^O^  ■^l'^  ^2 

-27i72cJ  +  ^  -27:72C 


2  2 


2  2 

7  +  7 

1  2 


(3.7) 


where  G  -  k2o^7i72.  If  1i~lz  then  G-g  defined  previously 


<7  -  y  <r"(ri)>  is  the  rms  height 


0  -  G(x,y)  is  the  surface  autocorrelation. 


•yrfy^yjT’jr 
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Substituting  (3.7)  in  (3.4) 


12 


X 


with 


r  -  G 


2  2 
7  +  7 

ll _ 2 


2^1 72J 


(3.8) 


The  integration  in  Rj^  can  be  performed.  From  integral  tables  (37,  p 
65. 6A) 


®  ik(R^-s)  -  ^(a2u2+/92V2) 

/  dR^De  -  I  e  ®  (3.9) 

>  00 


since  A  -  irUV 


Thus 


r 

1 2 


Ak  fi (e)f2 (e)e 


ik(r,  -r,  )_*r 


|^(a2U2+^2V2) 


8»r*r '  ^tj  Zi 


(3.10) 


2 .  Surface  Autocorrelation 

a .  Gaussian  Autocorrelation 


The  autocorrelation  in  this  case  is  known  to  be  (38) 


The  covariance  is  therefore 


r 

1  2 


-r-  ^  (a2U2+/32v2) 

A2k2f^  (e)f2 

Air^r'  r2 


•  I 

m— o 


m! 


Ni  ^ 


N, 


.  k2  ‘•'1  "2  2 


M1M2 


(3.15) 


We  shall  call  the  m-o  term  the  coherent  component.  It  is  equal  to 

-  ^(U'*(a2+4a2)+v2(^2+4^2)]+^l,(r  ) 

A2k2f,(e)f2(e)e  “  e'^ 

^12  Uit^v’  '^r^r2 

(3.16) 

The  same  result  is  obtained  using  (2.25)  and  (2.40),  calculating 

*  s,) 

<p(ri)>M<p(r2)>Me  - - - . 

If  the  surface  is  mirror- like,  G  approaches  zero.  Thus, 


r 

1 2 


M 


k2 


A^k^f^ (e)f2 (e)e 


[U2  (a2+4a2  ) +V2  (^2+4(02  ) 


ik(rj -r2 ) 


47r2r'  r2 


(3.17) 


Consequently,  (3.16)  can  be  re-expressed; 


r 

1 2 

c  o  h 


>  <p(ri)p*(r2)>j^e 


(3.18) 


so  that  the  degree  of  coherence  can  be  measured  by  e  ,  and  (3.15) 
is  expressed  in  terms  of  the  coherent  and  incoherent  components 


r  -r  e-" 

12  12M  m! 

m— 1 


-  k2  -  - 

with  a  factor  A  -  expi^Ca^U^+zaav®)  }  . 


.  2  rfN  1  2  fN  1  2 

_L  _2. 

m  ■  2  M,  M, 

o  uviy  vz</ 


(3.19) 


The  intensity  is  obtained  from  (3.19)  when  ri-r2 .  The  few 
changes  that  occur  in  writing  the  intensity  expression  are  worth 
enumerating,  since  they  will  be  invoked  a  few  times  throughout  thi 
s  tudy ; 

ri  -  rj  -  r 
71-72-7 

r  -  g 

a  “  0  a  —  -  € 

/3  -  0  ^  -  costf-5 

from  the  above,  fj^  (9)-f2  (9)  and  the  intensity  at  r 


.  ki  !i  '  ^ 

r  2M, 

UA  - — - — 

m-1 


with  N, 


and  r 


-fU  ^2  -  {cosd-S)V  A  -  exp(^(62u2  +  (cos^-5)2v2)  } 


A2k2  f  2  (9')  e 


-  ^[U2a2+V252 

a 

4w2(j-'r)2 


(3.20) 
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I  resembles  the  intensity  scattered  by  a  mirror- like  surface. 

1  1 M 

In  the  limit,  when  U  and  V  approach  zero,  the  expression  given  above 
equals  radiation  scattered  by  an  infinite  plane  surface  in  the  far 
field  --  see  Chapter  4.  The  first  term  in  (3.20)  is  the  coherent 
component  of  the  intensity: 


r  -re 

lleoh  IIM 


.0.21) 


2 .  Non-Gaussian  Surface  Autocorrelation 

a.  Slightly  Rough  Surface:  or  low  frequency  case, 
corresponds  to  GC«1.  The  exponent  involving  G  in  (3.10)  is  expanded 
in  a  series 


r 

12LF 


Ak2f  (0)f  (e)e^‘^^’^»''^2^e'^e 
1  2 

8w2r'2rir2 


^(a2U2+/92V2) 


X  /  drxDe^*'^^'®^l+GC+  ^  G^C^.  .  .  )  (3.22) 

The  first  term  is  familiar  to  us,  since  it  is  the  coherent  term 
that  appears  in  (3.16).  We  shall  be  interested  in  the  first  two 
terms  in  (3.22).  What  remains  therefore,  is  the  computation  of  the 
first  term  in  the  incoherent  part  of  the  covariance 


f  drxDe^’'^^'L'®^GC(rj.)  (3.23) 


!•  .  4 
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(3.26) 
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substituting  (3.29)  in  (3.10)  and  integrating  yields 


4ir2r'2ririBiB2 


-  |^[U2 (02+4^2/62 )+v2 (^2+4^ /B^)]  -  |ro 
X  e  e 


(3.30) 


for  covariance , 


where  B^  -  yi-GTiU^ 


B,  -  yi-GT,V2. 


Note  that  and  T2  are  negative,  B^  and  B2  are  positive  roots. 


For  intensity  (3.30)  gives 


A2k2f2(e)e 


-  ^[(?/Bf)U2+(^/B^)V2l 


4,r2(r'r)2B2B2 


(3.31) 


B.  THE  SLOPE -OPERATOR  MODEL 


1 .  Formulation 


This  model  shall  be  developed  using  both  the  Fresnel  and 


Frauhofer  approximation.  Firstly,  u.sing  (2.36)  the  pressure  at  some 


point  r2  is  derived 


-  2.ir'r, 


(*2^'  '  '  '  ’  -sin^2  )dx'  ’  ’  dy  '  ' 


(3.32) 


Multiplying  (2.36)  by  the  complex  conjugate  of  (3.32)  yields 


52 


j^2gik(ri  -r_) 


P(ri)P*(r2)-  /  ///;  •■*"') 

12 


■‘■^25‘y  -sintf2  ]dx- ■  dy  -  dx’ •  ■  dy 


(3.33) 


and  the  expectation  of  (3,33),  becomes 


,  ,  ik(r,  -r,  )  ®  . 

<p(r,  )p.(r,  »  -  SffI  D.  .  D.  • .  ■•»'•> 


■‘■^afy  -sintfj  ]>dx' '  dy  ■  dx- ■  ■  dy  '  ' 

(3.34) 

The  random  quantity  f  is  in  the  phase,  therefore  in  the  center-of- 
mass  and  difference  coordinates 

<eik(*'  ■-*"')>  _  ^ik(rj_-s+Rj.-s]^^-ik(7ir'  '  -Tzf  '  ’  )> 


x^+v^  -  - 

^  •  • 

p 

4  +  x2+y2 

+  P  (xy+xy) 

(3.35) 


where  -  P1+P2  P'  -  P1-P2 

The  first  two  exponentials  are  the  familiar  Fraunhofer 
approximation  and  characteristic  function.  The  third  term  is  the 
correction  necessary  so  that  the  Fresnel  phase  approximation  can  be 
presented.  This  term  is  the  first  in  the  expansion  to  show  coupling 
among  the  orthogonal  coordinates . 


Welton  (25)  and  others  (39,40)  apply  a  theorem  of  mathematical 


statistics  (41,  p.  68)  that  postulates 


3 _  .  -ikbz^  3  -ikbz^ 

~ —  <e  >  -  e  > 

3q  3q 


(3.36) 


where  z-z(q),  provided  the  first-  and  second-order  derivative  of  the 
autocovariance  exist. 


Thus , 


__L  S — -72f  ’  "  ) 
lk7i  3x"  ® 


■■  ■••^-ik(7if-72r-")>  .  ^ _i —  <g-iK(7ir'--72r'’')> 

'x  ^*71 72  '  Sx"  ' 

(3.37) 


Application  of  this  theorem  allows  expressing  the  averaged  quantity 


in  (3.34)  as 


ik[p-(^^^^  +  x2+y2)+p+(xy+xy)]  r  “  o  i  -w  ,  ^  n 


(3.38) 


where  the  Slope-Operator  is  defined  as 


T( 


X.  ■  ,x-  •  •  .y  •  .y  ■  • 


■•■^1^2  Qy,  ,  3y.  ,,"*"*i^2  3x' '  ay  '  ■  ^2  Sx' '  ’  ay  ' 


ik  [7isin«,[  52f^,]-72sinS2[  6,  ]]} 


(3.39) 


in  the  old  coordinates.  In  the  new  coordinates 


T(x.y)-sintfisin«2- 


SI- 


*  (  £ 


^  ^  .32 _ 

l*2+«2*l)axay 


+ik 


:  7iSintfi  ^ 


£  ^ 
2  ax 


«2fy)-hr2sintf 


'  f  £  ^ 

2[  ‘lax 


-4)  ]} 


(3.41 


Incorporating  all  these  changes,  as  well  as  replacing  the 
characteristic  by  eq.  (3.7),  transforms  (3.34)  into 

P’  (^  +x^+y2 )+P'''(xy+xy ) 


r  - 

12 


2  ik(ri-r2)  -r  « 


k^e 


X  iVc 


4x2r ' ^r, r. 


J  drj^  f  dRj^  DDe 


^  gik(rj^s+Rj_-s).^^GC(r^) 


(3.41) 


The  first  exponential  within  the  integral  is  identically  unity 


when  the  Fraunhofer  phase  approximation  is  discussed.  Since  T  is  an 
operator  in  rj^  then  the  integration  with  respect  to  Rj^  can  be 
performed  using  the  result  from  (3.9) 


r 

12 


FH 


Ak£e 


8tr=*r'27i7, 


r^j  IN  ikr  I  •  s_  GCCr, ) 
I  dr^  De  Te 


(3.42) 

for  the  Fraunhofer  approximation  case.  Henceforth  the  superscript 
Fh  indicates  Fraunhofer  and  the  superscript  Fn  indicates  the  Fresnel 
result . 

The  Fresnel  result  is 


Fn 


kUs  2  ^  klJi£^  2  = 


12 


2jr2r'2rir2Wi 


«  ik[P'(  )+P'*'xy] 

X  J  drj^De  e 


4  -^Vikr^.s^^GC(r^) 


(3.43) 


with  Qi  -  y2-ikP*U2  Q2  -  JZ-ikF'V^ 

-  74Q^Q^  +  (kUVP‘'') *  ,  all  positive  roots 

3 .  Surface  Autocorrelation 

3 . 1  Surface  with  Gaussian  Autocorrelation 

The  autocorrelation  for  this  case  appears  in  (3.11)  and  the 


result  of  the  slope  operator  acting  on  (3.12)  is 


.-  'jr"j."»,  r". liTV’*  y 


j  u  i  -P  ak2  - 

and  b  -  -  -ik^) 


The  coherent  component  in  either  approximation  is 


r  -  r  e' 

12  coh  12  M 


j  —  Fh  or  Fn 


(3.47) 


and  the  covariance  from  a  mirror- like  surface  is 


-  '^[U2(a2+4a2)+V2(^2+4^2)| 


4jr2rir2r'2 


s in^  j  s 10^2 


(3.48) 


for  the  Fraunhofer  case,  obtained  by  taking  the  m-0  term  from  (3.45) 


and  letting  a  and  L^,L2  approach  zero  and  infinity,  respectively. 


For  the  Fresnel  approximation  we  use  (3.46)  and  perform  the  identical 


procedure 
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To  get  the  intensity  we  let  equal  to  r2 ,  (3.45)  and  (3.46)  for  the 
Fraunhofer  case, 
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Again,  the  coherent  component  of  the  intensity  is 
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the  covariance  can  be  expressed  in  terms  of  its  coherent  and 
incoherent  component  as 
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where  the  argument  in  square  brackets  [  ]  is  that  of  (3.45). 
Similarly,  the  intensity  is 
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The  first  term  in  both  expressions  is  their  respective  coherent 
component . 

b .  Non-Gaussian  Surface  Autocorrelation  in  the  Farfield 
1 .  Preliminaries 

The  Fresnel  phase  approximation  result  shall  not  be  considered 
herein,  hence  there  is  a  need  to  specify  that  source  and  receivers 
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be  in  the  farfield  and  the  beawidth  be  sufficiently  narrow.  The 
criteria  for  farfield  for  the  Tangent  Plane  methods  was  developed  and 
discussed  by  Welton  (25)  and  can  be  summarized,  for  the  worst  case 
(specular  direction)  intensity  computations  he  finds  that 


2  >  UM  (1,  lEl, 

L  ~  -1/2  ri  '  sinBi 


(3.57) 


In  the  nonspecular  direction  the  criteria  is  less  stringent. 

For  covariance,  both  receivers  must  satisfy  the  above  criterion,  thus 
r  and  0^  must  assume  the  values  of  the  respective  receiver. 
Furthermore,  we  limit  ourselves  to  presenting  results  valid  for  Che 
surfaces  with  onmidirecCional  spectra,  hence  Li-L2-2‘ ^ '  ^L. 

The  same  high  and  low  frequency  approximation  performed  in 
connection  with  the  Standard  Model  shall  also  be  used  here.  The  term 
involving  the  operator  and  the  surface  autocorrelation  needs  to  be 
computed  for  both  the  high  and  low  frequency  approximations : 
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where  S(  ) 


-/  dxW(x)(  )e^^‘''^ 


and  the  quantity  in  parenthesis  is  always  a  scalar  (magnitude) . 
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2.  Slightly  Rough  Surface:  The  exponential  in  (3.58)  is 
expanded.  The  first  2  terms  are  retained  and  (3.58)  is  substituted 
into  (3.42),  the  equation  corresponding  to  the  Fraunhofer 
approximation  result.  The  integrations  in  rj^  are  peformed  to  obtain 
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The  range  of  validity  of  (3.59)  is  that  GC(rj^)«l.  The  intensity  is 


r  -  /  2/  ,  {k272sin2^.e 

11  LF  4jT2(r'r)*72  ^ 


-  |^[a2U2+^2v2' 


+  2^02  i"  dxW(x)e 


(ka+x^)^  (kp+xy)^ 
- —  U2-  V2 


X  (k272sin2tf^+e2j^2^52^2^2€5x'xXy+2k7sin«i  ( ejfx+'^X’y) ) 


+  /  dX  ;  dx'  W(x)W(x')e 


( ka+Xx+X^ )  ^  ( k/3+Xy+x„ )  ^ 

— ^ ^  U2-  - V2 


(X^+XxX^)+« ^  (Xy+XyXy)+2€«Xx(Xy+Xy) 
+2k7sin5i (exx+^Xy)) 


CO  00  CO 


+  (2xg'2)3  /  dx  /  dx'  J  dx"  W(x)W(x')W(x") 


•  CO  -  CO  >00 


( ka+Xx+X„+X.j, )  ^  ( k^+Xy+X.^,+X.^^ )  ^ 

- ii - a  It2.  -^- - ^ i  V2 

2  2 


(f ^XxX^+'J^XyX  +2e5xxX  ) )  (3.60) 


3.  Very  Rough  Surface:  The  method  of  stationary  phase  is 
employed  to  evaluate  (3.42)  with  respect  to  its  spatial  dependence 
once  (3.58)  is  substituted  in  (3.42)  and  the  results  from  (3.28)  and 
(3.29)  are  used.  In  this  case,  the  result 
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where  -  yi-GTj^U^  and  B2  ”  J\-GT2^^  sre  positive  roots. 
This  formulation  is  valid  when  GC»1  and  the  intensity 


r 

1 1 


Fh 

HF 


A^e  ^ _ 

BiB24x2(r'r)272 


{k^y^sin^s^e 


kifalUi 

2  B2 


B2 

2 


+ 


I 


d;t«(Ar)e 


(ka+Xx)^U2 

2B2 

1 


(k^+Xy)2V2 

2B2 

2 


66 


(f  ^Xx'*'^^Xy+2«5)fxXy-*-2k7sin^i  (eXx+^X’y)  '> 
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Low  and  high  frequency  expressions  in  the  Fresnel  approximation  could 
be  developed.  The  algebra  is  involved  and  tedius ,  but  not 
impossible . 

C.  THE  COMPOSITE -ROUGHNESS  MODEL  (5,6,7) 

1 .  Introduction 

The  present  section  deals  with  a  technique  that  has  been 
successful  in  studies  involving  a  one -dimensional  ocean  surface. 

This  is  so  because  it  assumes  that  the  surface  to  which  this  model 
applies  scatters  the  incoming  wave  in  two  phenomenologically 
different  ways.  First,  the  long  wavelength  incoming  radiation,  sees 
a  surface  that  is  locally  flat,  possibly  with  non-zero  slope,  and 
scatters  from  it  as  if  though  it  were  a  distribution  of  mirror- like 
facets.  In  optics  this  corresponds  to  the  physical  optics  domain  of 
radiation  scatter  and  is  not  new  to  us  since  it  is  formulated  by  the 
Helmholtz -Kirchhoff  Integral  or  the  Double-Layer  Potential  Integral 
discussed  previously.  Second,  the  short  wavelength  incident  wave 
will  be  scattered  diffusely  by  the  ripple  which  rides  on  the  large 


scale  undulating  surface  (see  Fig.  3).  The  ripple  has  the  features 
of  surfaces  to  which  the  Method  of  Small  Perturbation  applies ,  such 
as  rapidly  changing  slopes  and  small  crest- to- trough  amplitudes.  A 
theory  that  applies  to  surfaces  of  this  type,  with  two  roughness 
ranges,  could  be  one  that  merges  the  Helmholtz -Kirchhoff  Integral  and 
the  Method  of  Small  Perturbation.  It  would  imply, that  the  scattered 
wave  undergoes  Dali-like  distortion  as  predicted  by  the  Helmholtz 
Integral,  but  the  sharp  definition  is  somewhat  lost  due  to  the 
contribution  of  the  perturbation. 

The  model  being  discussed  is  the  Composite -Roughness  Model.  It 
applies  well  to  surfaces  that  can  be  divided  into  two  scales  of 
roughness,  and  have  most  of  their  energy  scattering  potential 
concentrated  in  the  lower  regime  of  the  spectrum.  However,  a  problem 
arises  when  it  is  not  clear  where  this  spectrum  splitting  occurs. 
McDaniel  (7)  has  improved  the  method  by  including  diffraction 
corrections.  These  corrections  are  also  beneficial  in  reducing  the 
sensitivity  of  the  results  to  changes  in  cutoff  surface  wavenumber. 

As  with  the  Tangent-Plane  methods  this  strategy  requires  the 
computation  of  the  expectation  of  the  projection  of  the  source 
wavenumber  and  the  random  surface  normal.  A  choice  must  be  made, 
which  usually  considers  computational  simplicity  heavily  on  the 
averaging  technique.  But  the  demands  here  are  not  so  stringent  as  in 
the  Tangent- Plane  Method  since  the  normal  in  question  is  taken  to  be 
the  slowly  varying  large-scale  roughness  normal.  The  next  two 
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subsections  will  present  the  spatial  covar-.ance  formulation  using  the 
conventional  Composite-Roughness  Model  and  diffractive  corrections 
are  introduced  in  the  last  section. 

2.  The  Conventional  Model:  the  interest  is  to  develop  two 
expressions  in  the  Fraunhofer  phase  approximation.  One  will 
incorporate  the  function  f(0)  as  an  approximation  to  the  evaluation 
of  the  dot  product  involving  the  random  normal  to  the  surface. 

Another  expression  will  use  the  slope  operator  instead.  The  last 
section  shall  leave  the  issue  of  the  dot  product  behind  and  shall 
describe  the  changes  that  occur  in  the  model  when  diffractive 
corrections  as  proposed  by  McDaniel  (7)  are  incorporated. 

a.  Inteeration-bv-Parts  Procedure:  Equation  (3.10)  is  our 

starting  point.  The  autocorrelation  function  is  expressed  as  a  sum 
of  contributions  from  the  spectra  of  the  large  scale  surface  plus  the 
small  scale  surface  spectra. 

C(rj^)  -  CL(r^)  +  Cs(rj^)  (3.63) 
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expression  for  the  covariance  is 
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The  large-scale  surface  autocrrrelation  can  be  expanded  about  rj^-0 . 
The  autocorrelation  length  is  assigned  to  be  small  and  the  spectrum  is 
assumed  non-directional .  Thus, 


(3.66) 


with  CTL  ■  rms  large-scale  surface  waveheight 
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■  the  X  component  of  the  rms  slope -large  scale  surface 
«  the  y  component  of  the  rms  slope -large  scale  surface. 

Substitution  into  the  explicit  form  of  (3.65)  gives 
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The  symbol  |  indicates  that  the  integrals  lower  limit  corresponds  to 
the  cutoff  spectral  wavenumber. 

Integrating  over  x  and  y  gives 
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are  positive  roots. 


and  the  intensity  is  similar  except  that  the  constant  in  front  of 
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resultant  equation  for  intensity  are  valid  so  long  as  (~^)^  «1. 

In  terms  of  a  coherent  and  incoherent  component,  (3.68)  can  be 
expressed  as 
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The  intensity  has  the  same  term  in  (  )  brackets  (the  values  of  a  and 
/9  are  different  for  intensity)  but 
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b .  Incorporating  the  Slone  Operator: 

Recognizing  the  similarities  between  the  expansion  in  (3.64)  and 
the  one  used  to  obtain  (3.59),  the  result  of  using  the  slope -operator 
in  the  context  of  the  Composite -Roughness  Model  can  readily  be 
derived  using  the  results  of  (3.68): 
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valid  so  long  (^^)«1.  The  symbol  ^  indicates  that  the  lower  limit 


of  integration  is  consistent  with  the  lower  limit  in  (3.68). 
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The  surface  correlation  function  is  partitioned  as  was  done  in 
(3.63),  but  the  expression  for  the  large-scale  surface  correlation  is 


..  .  Xj^  2x 

ct^Cl  ”  *  2  ^  /  dnW(x,fl){cos(xTj^) 

^  o  o 


•1+  i(x-ri)2} 


(3.74) 


Using  (3.63)  and  (3.77),  the  exponential  term  in  (3.10)  is 

—  r*  f  f  y  n 

ik(rj_-S)-  [(a  v)2+(o  ^y)2]+-^  /KcU/d£lW(x .Q)  {cos(x- r^) -1 

2<7*  ^  2xa*  o  o 


(3.75) 


When  G(~)  is  small  the  above  result  can  be  recast  as 

<j 


{ik(rj^-s)-  [(<7j^x)2+(a^)2])[l-G2(l-C(rj^))+  ( (c^^^x)  *+(ct^)  * )  ]  . 

(3.76) 

Using  (3.75),  or  (3.76)  when  applicable,  in  (3.10)  and  performing  the 
integrations  over  rj^  variables  yields  an  expression  that  incorporates 
diffractive  corrections  and  that  is  less  sensitive  to  the  choice  of 
partitioning  wavenumber.  The  result  is  also  attractive  because  it 
lacks  computational  complexity  and  gives  fairly  good  agreement  when 
its  predictions  are  compared  to  experiment. 


CHAPTER  IV 


SCATTERING  STRENGTH  AND  COHERENCE 


A.  SCATTERING  STRENGTH 

The  dlfferenc  models  presented  in  the  previous  chapter  are 
compared  and  their  behavior  analyzed.  The  aim  of  the  first  half  of 
this  chapter  Is  to  compare  the  Standard  and  Slope -Operator  Model  In 
general  terms,  In  addition  to  the  three  alternatives  to  the  Slope- 
Operator  Model:  the  Fraunhofer  Approximation,  the  Fresnel 

a 

Approximation  and  the  case  when  the  random  normal  ^  Is  approximated 
by  ^  .  To  enable  comparison,  section  A  examines  the  scattering 
strength,  defined  as 


s  -  ioiog„ 


(^.1) 


Involving  a  flat  and  a  Gaussian  surface.  Note,  however,  that  the 
plots  shown  are  normalized  to  the  maximum  values.  Therefore,  the 
figures  show  the  normalized  scattering  strength.  The  second  half  of 
the  chapter  compares  the  Standard  and  Slope-Operator  models,  as  they 
are  applied  to  the  measure  of  coherence  (1,  p.  498) 
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when  the  surface  Involved  is  Gaussian.  From  the  three  alternatives 


of  the  Slope-Operator  Model,  the  Fraunhofer  Slope -Operator  Model  is 
chosen  for  comparison  to  the  Standard  Model. 


The  classical  expression  for  the  scattering  strength  Involving  a 
flat  surface  In  the  Fraunhofer  approximation  Is 


Ak^sin^tf 

So  -  ‘l 


(4.3) 


and  is  shown  in  Figure  4  in  the  specular  case.  Equation  (3.51),  the 


Slope-Operator  Model  yields 

Ak^sin^tf 

S  -  lOlog  {  -r-^ - ^). 

•  o  *10  4ir2 


(4.4) 


when  r'  and  r  become  large  and  the  beam  very  narrow.  If  the  same 
approximations  were  possible  for  the  Standard  Model  (3.20)  the  result 
would  be 


s.  -  lOlog,, 


(4.5) 


For  the  specular  test  shown  in  Figure  4,  f(0)  reduces  to 
(l+sin^tf )/2sinfl .  Figure  5  shows  good  agreement  between  the 
Fraunhofer  Slope -Operator,  Standard  and  Figure  4.  The  plots  of 
Figure  5  show  the  models  without  making  any  approximations.  The 
Fresnel  expression  will  agree  with  the  other  two  if  the  above 
mentioned  limiting  vaues  for  r' ,  r  and  the  beam  are  made.  The 
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Figure  4.  Scattered  intensity  specular  direction  in  the  Fraunhofer 
approximation.  Flat  surface.  Source  at  500  cm  frequency 
is  10  KHz ,  beamwidth  is  3" . 


discrepancy  between  the  Standard  Model  and  the  Slope -Operator  Model 
in  connection  to  the  angular  dependence  demonstrates  how  weakly  they 
depend  on  the  functions  f(0)  and  sin^0  in  the  specular  case. 


For  the  forward  case,  the  lack  of  consistency  between  the  two 
models  becomes  more  obvious.  This  hypothesis  is  based  not  only  on 
the  fact  that  f(8)  does  not  simplify  to  sln^^j^  for  all  cases,  but 
because  by  assumption  the  Standard  Model  is  valid  not  far  from  the 
specular  direction.  Figure  6  illustrates  this  point.  The  Standard 
Model  misbehaves  at  low  and  very  high  receiver  grazing  angles. 

Since  the  Standard  Model  is  limited  to  the  near  specular 
direction,  in  the  backscatter  test  it  is  expected  to  perform  poorly 
far  away  from  the  normal  to  the  mean  surface.  This  is  obviously 
demonstrated  in  Figure  7 .  These  figures  show  the  backscatter 
involving  a  flat  surface  and  a  slightly  rough  Gaussian  surface  for  a 
source  with  3°  beamwldth  operating  at  10  KHz.  In  this  test  it  is 
expected  that  the  energy  be  higher  in  amplitude  and  largely 
concentrated  around  the  normal  (source  grazing  angle  of  90®)  for  the 
flat  surface  scatter  as  shown  in  Figures  7a  an  7b  for  the  Fresnel  and 
Fraunhofer  Slope -Operator  Models  respectively.  This  is  not  the  case 
for  the  Standard  Model  thac  is  plotted  in  7c.  This  last  figure  shows 
good  agreement  in  the  normal  direction  but  breaks  down  for  low 
grazing  angles,  because  the  model  is  used  outside  its  region  of 
validity  (see  II.B.l). 
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2 .  Gaussian  Surface: 

a.  Backscatter 

Figures  8a,  8b,  and  8c  correspond  to  equations  (3.51),  (3.50), 
and  (3.20),  respectively,  in  the  backscatter  case.  The  three  figures 
show  good  agreement  in  the  prediction  of  backscatter  close  to  the 
normal  direction.  As  expected,  the  low  grazing  angle  regime  shows 
the  Fresnel  expression  yielding  demonstrably  less  backscatter.  Near 
the  normal  the  three  models  show  that  the  amplitude  changes  in 
proportion  to  the  frequency  squared,  which  is  implicit  in  the  k^  that 
is  part  of  the  coherent  component  of  equations  (3.51),  (3.50),  and 
(3.20) . 

The  effect  of  changing  the  rms  height  of  the  surface  is 
illustrated  in  Figures  9a,  9b,  and  9c.  With  an  increase  in  the  rms 
height,  the  area  under  the  curve  becomes  wider,  as  is  the  case  in  9a 
and  9c,  and  the  overall  amplitude  drops  as  seen  in  all  three  plots. 
Again,  the  Standard  Model  fails  to  provide  good  results  far  from  the 
normal  direction. 

Further  insight  into  the  way  the  Standard  Model  behaves  for  low 
source  grazing  angles  is  given  by  Figures  10a,  b,  c.  In  these  three 
figures  the  broken  lines  represent  the  coherent  component  of  the 
backscattered  intensity  and  they  can  be  compared  to  the  solid  line 
which  represents  the  coherent  and  incoherent  components  combined.  In 
the  Slope-Operator  Models  the  incoherent  component  contributes  to  the 
radiation  even  at  low  grazing  angles.  This  is  not  the  case  in  the 


Frequency  dependence  of  backscattered  radiation 
involving  Gaussian  surface.  (a)  Fresnel  slope- 
operator,  (b)  Fraunhofer  slope-operator,  (c)  standard. 
_ :  5KHz,  - :  7KHz.  ****:  12KHz 


Figure  8 . 
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Figure  9. 


Backscatter.  Variation  with  nns  height. 

_ ;  0.5cm,  - :  1.0cm,  ****:  2.0cm. 

(a)  Fresnel  slope-operator,  (b)  standard, 
(c)  Fraunhofer  slope-operator 


Standard  Model,  for  which  the  radiation  is  exc] isively  predicted  by 
the  coherent  term.  The  plots  correspond  to  a  source  operating  with  a 
3“  beamwidth  at  a  frequency  of  10  KHz.  In  all  the  models  the 
incoherent  component  is  the  stun  of  ten  terms  which  is  more  than 
adequate  for  the  values  to  converge  for  this  frequency.  In  general, 
the  higher  the  frequency,  the  more  incoherent  texrms  are  needed  since 
an  expansion  of  the  exponential  involving  the  roughness  parameter  G 
was  performed  in  all  the  models.  The  result  in  turn  depends  on  the 
square  of  the  frequency. 

b.  Forward  Scatter 

In  the  forward  scatter  test,  a  source  is  placed  at  500  cm  from 
the  surface,  at  an  angle  of  45“ .  The  receiver  is  200  cm  away  from 
the  surface.  Again,  the  surface  has  autocorrelation  lengths  of  20  cm 
and  an  rms  height  of  2  cm.  The  scattering  strength  is  expected  to 
peak  in  the  neighborhood  of  a  receiver  grazing  angle  of  45“ , 
corresponding  to  the  specular  position. 

Figure  11  is  the  forward  scatter,  for  a  source  with  30“  beam 
operating  at  10  KHz,  as  predicted  by  the  three  models.  The  peak 
occurs  in  the  neighborhood  of  48“ ,  not  at  the  specular  angle.  The 
peak  is  frequency  and  beamwidth  dependent  in  all  three  models  and 
will  shift  to  the  specular  angle  as  the  frequency  or  the  beamwidth  is 
increased.  This  shift  is  apparent  in  Figures  12a,  12b,  and  12c,  in 
which  all  models  exhibit  a  peak  which  is  further  away  from  the 


specular  angle  at  lower  frequencies. 
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Figure  11.  Forward  scatter  for  a  source  at  500  cm  and  45°  with  a 
3°  beam  and  an  operating  frequency  of  10  KHz . 

_ :  Fresnel  slope -operator,  - :  standard, 

****:  Fraunhofer  slope -operator 


Figures  13a,  13b,  and  13c,  show  the  maxims  moving  away  from  the 
specular  angle  for  smaller  beamwidths.  The  plots  shown  in  Figures 
14a,  b,  c  hint  that  the  problem  lies  in  the  coherent  component  of 
equations  (3.20),  (3.56),  and  (3.51).  Since  it  involves  frequency 
and  beam  characteristics  and  is  a  function  of  geometry,  the 
exponential  in  the  coherent  term  must  be  causing  the  phenomenon. 

The  coherent  components  of  all  three  models  are  compared  in 
Figure  15  to  show  once  more  that  the  Standard  Model  misbehaves  for 
angles  far  away  from  the  specular  direction  (80°).  This  is 
particularly  obvious  for  shallow  receiver  angles.  A  comparison  of 
this  figure  with  number  11  gives  an  idea  of  the  importance  of  the 


incoherent  component  away  from  the  specular  direction. 


In  this  test  the  source  and  surface  are  as  in  the  previous 
section,  and  the  receiver  is  at  200  cm  from  the  surface,  the 
azimuthal  angle  is  allowed  to  change  while  the  depression  angle  is 
held  fixed  at  45° ,  The  reader  is  reminded  that  the  forward  direction 
corresponds  to  an  azimuthal  angle  of  90° ,  whereas  the  side  and  back 
are  at  180°  and  270°  respectively,  as  depicted  in  Figure  2.  Figure 
16  shows  a  comparison  of  the  three  models.  The  peak  occurs  at  90° 
which  corresponds  to  the  specular  direction,  facing  the  source.  The 
dip  at  270° ,  in  turn,  is  the  antispecular  or  source  position.  The 
three  models  tend  to  behave  similarly,  although  the  amplitudes  do  not 
agree  except  for  when  the  receiver  is  in  the  specular  direction  where 
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Figure  13. 


Forward  scatter,  (a)  Fresnel  slope - operator , 

(b)  standard,  (c)  Fraunhofer  slope -operator . 

.  50^  _ .  go ^  ***★;  10°  beamwidth  for  a  source 

operating  at  lOKHz 
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Azimuthal  forward  scatter.  Source  at  500  cm,  45® 

grazing  angle.  _ :  Fraunhofer  slope-operator, 

standard,  ****:  Fresnel  slope-operator. 


Figure  16. 


the  Standard  and  Fresnel  result  show  agreement.  Figures  17a,  b,  c 
show  general  qualitative  agreement.  Changes  in  the  symmetry  of  the 
surface  affect  the  predicted  scatter  in  the  antispecular  direction 
(270").  When  the  autocorrelation  length  parallel  to  the  surface 
direction  is  increased,  the  backscatter  drops  drastically  towards  the 
source  and  the  forward  scatter  increases  somewhat.  In  the  limit  of 
large  autocorrelation,  this  result  is  logical  since  in  the  mirror 
surface,  the  scatter  is  mainly  in  the  forward  direction.  On  the 
other  hand,  when  the  cross-directed  autocorrelation  length  is 
increased  the  figures  display  a  narrowing  of  the  forward- directed 
scatter  and  an  increase  in  scatter  towards  the  source.  If  a  surface 
with  one -dimensional  corrugation  is  envisioned  and  a  source  is 
directed  along  the  corrugation  we  would  have  the  limit  of  the  first 
example.  When  the  source  is  directed  perpendicular  to  the 
corrugation  we  would  have  the  limit  of  the  second  example. 

d.  Summary 

The  Standard  and  Slope -Operator  models  are  qualitatively  alike, 
except  for  low  grazing  angles.  The  energy  distribution  is  better 
described  by  the  Slope-Operator  Model  than  the  Standard  Model  as  the 
surface  parameters  change,  mainly  because  the  energy  in  the  Standard 
Model  is  exclusively  dominated  by  the  coherent  contribution.  The 
coherent  component  is  a  good  roughness  estimator.  A  different 
approach  to  the  integration-by-parts  or  the  Slope -Operator  is 
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Figure  17. 
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Azimuthal  forward  scatter.  :  20,20,  40,20, 

*★**:  20,40cra  for  the  x  and  y  directed  autocorrelation 
lengths  respectively.  (a)  Fraunhofer  slope -operator , 
(b)  standard,  (c)  Fresnel  slope-operator. 
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possible  when  the  surface  is  only  slightly  rough.  Th^s 
is  done  by  approximating  ^  upon  examining  the 

development  in  analogy  to  the  Slope-Operator  Model,  assuming  and 
fy  to  be  zero  in  eq.  3.34.  The  result  is  identical  to  the  coherent 
component  of  3.45.  Thus,  a  comparison  between  any  of  the  previously 
mentioned  models  with  this  new  approach  is  one  in  which  the  coherent 
component  of  the  Slope -Operator  Model  is  compared  with  the  full  model 
as  well  as  the  Standard  approach.  The  main  deficiency  of  the  latter 
model  is  that  no  spectral  information  is  supplied,  or  required,  by 
the  resulting  expression.  Hence,  the  approximation  is  valid  in 
regions  where  the  field  is  coherent.  That  comprises  a  region  in  the 
neighborhood  of  the  specular  direction  for  forward  scatter  and  in  the 
normal  direction  for  backscatter.  The  size  of  these  regions  grows  as 
the  autocorrelation  lengths  get  large  and  the  rms  height  gets  small. 

The  significant  difference  between  the  two  Slope-Operator  Models 
in  the  far- field  occurs  at  low  grazing  angle  where  the  Fresnel 
expression  predicts  lower  energy  levels.  To  answer  the  question  of 
which  model  will  agree  better  with  experiment,  a  comparison  with  data 
is  required.  In  the  following  section,  among  the  Slope -Operator 
formulations,  the  Fraunhofer  expression  is  chosen  to  study  coherence. 
B .  COHERENCE 

A  source  operating  at  10  KHz  with  a  3°  beam  is  placed  at  200  cm 
and  45°  with  respect  to  the  mean  insonified  surface.  The  surface 


has  autocorrelation  lengths  of  20  cm  and  an  rms  height  of  2  cm.  When 
one  receiver  is  placed  at  the  source  position  and  another  is  allowed 


to  move  polarly  as  in  Figure  18a,  the  coherence  magnitude  varies  as 
shown  in  Figure  19a  when  predicted  by  the  Standard  and  Slope-Operator 
Models.  The  peak  occurs  at  135°  which  corresponds  to  the 
source/fixed- receiver  location.  If  the  moving  receiver  travels 
azimuthally  at  a  constant  depression  angle  of  45°  as  per  Figure  18b, 
the  magnitude  coherence  is  predicted  as  plotted  in  19b.  The 
coherence  is  unity  at  the  source  location,  when  the  azimuthal  angle 
is  270°. 

When  the  fixed  probe  is  placed  200  cm  away  from  the  surface  with 
a  depression  angle  of  45  and  a  moving  receiver  travels  vertically 
forming  a  plane  perpendicular  to  the  mean  surface,  the  arrangement  is 
called  here  a  vertical  coherence  test.  Figure  20a  illustrates  this 
arrangement.  The  angle  ^  is  set  at  90°  for  the  forward  scatter 
direction,  180°  for  the  side  scatter  direction  or  270°  for 
backscatter  direction.  The  results  for  these  three  cases  appear  in 
Figure  21a,  b,  c  respectively.  The  peaks  occur  at  a  horizontal 
distance  of  141.42  cm,  which  corresponds  to  the  depth  of  the  fixed 
probe  (i.e.,  141.42-sec  45°-200.0).  Since  r^  is  different  from  r2 
except  at  the  probe  location,  the  phase  of  the  coherence  will  be  non¬ 
zero  as  shown  in  Figure  22.  This  last  figure  simply  reflects  the 
change  in  phase  due  to  the  k(r^-r2)  imaginary  argument  in  the 
exponential  of  eq.  3.45  for  the  Slope-Operator,  and  3.19  for  the 
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Vertical  coherence.  _ ;  standard,  :  slope- 

operator.  (a)  forward  scatter  direction,  (b)  side 
scatter  direction,  (c)  back  scatter  direction. 
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Standard  Case.  The  horizontal  coherence  test  is  illustrated  by 
Figure  20b.  The  fixed  probe  is  placed  at  a  45“  depression  angle,  200 
cm  away  from  the  surface.  The  forward,  side,  and  back  directions 
correspond  to  azimuthal  angles  90“ ,  180“ ,  and  270“  respectively  (see 
Figure  2) .  The  moveable  receiver  travels  in  straight  line  parallel 
to  the  mean  surface  within  the  azimuthal  angles  shown  in  Figures  23a, 
b,  c.  The  coherence  is  strongly  dependent  on  exp{-ri  from  expression 
(3.45)  and  (3.19)  which  is  a  function  of  source  and  receiver 
depression  angles  as  well  as  surface  roughness.  These  are  factors 
which  strongly  affect  the  vertical  coherence,  but  minimally  affect 
the  horizontal  coherence,  explaining  why  so  little  information  is 
provided  by  the  horizontal  measurement  while  the  vertical  coherence 
is  a  good  indicator  of  source  direction  as  well  as  surface  roughness. 

The  width  of  the  magnitude  of  the  coherence  tends  to  get 
narrower  with  an  increase  in  frequency  as  shown  in  Figures  24  a,b. 

The  figures  show  that  both  models  tend  to  better  agreement  as  the 
source  frequency  is  increased.  On  the  other  hand,  better  agreement 
between  both  models  occurs  when  the  beamwidth  is  decreased.  Figures 
25  a,b  are  plots  of  the  coherence  as  a  function  of  beamwidths  for  the 
two  models.  The  coherence  curves  become  narrow  with  an  increase  of 
frequency  or  beamwidth  due  to  the  exponentials  involving  -k^  as  well 
as  -U^  and  -V^ . 
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Figure  23. 
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Horizontal  coherence  _ :  standard,  - :  slope- 

operator  (a)  backscatter  direction,  (b)  side  scatter 
direction,  (c)  forward  scatter  direction.  Note 
extremely  high  resolution  of  vertical  axis  when 
comparing  these  plots  to  vertical  coherence  plots. 
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Coherence  frequency  dependence.  _ :  5KHz ,  - ; 

lOKHz ,  ****:  12KH2  (a)  standard  model,  (b)  slope- 
operacor  model 


Figure  24. 
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The  cc-herence  is  a  function  of  source  position.  For  example, 
when  the  source  grazing  angle  is  shallow,  is  large  so  that  the 


coherence  is  low.  This  is  particularly  tue  when  a  and  a  are  non¬ 
zero  as  is  the  case  when  the  source  and  receivers  are  non-coplanar . 
Coherence  is  sensitive  to  variation  of  the  rms  height  but  is  best 
used  as  an  estimator  of  spectral  information.  The  variation  of  the 
coherence  with  rms  height  is  shown  in  Figures  26  a,b.  Figures  26  c,d 
show  that  increasing  the  autocorrelation  length  parallel  to  the 
source  direction  will  tend  to  broaden  the  coherence  curves.  Changing 
the  cross -directed  autocorrelation  length  does  not  have  any  affect  on 
the  coherence  plots,  as  is  evident  in  Figures  27  a,b.  The  best 
estimator  of  surface  roughness  is  the  coherent^component  of  the 
intensity,  since  it  is  equal  to  <p>„<p>*^e  ^  ,  i.e.,  the 

coherent  intensity  reflected  from  a  mirror  surface  times  a 
"coherence"  factor  which  is  a  function  of  surface  roughness. 
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Figure  26 . 


Coherence.  (a)  dependence  on  surface  inns  height, 
standard  model,  (b)  dependence  on  surface  rms  height, 

slope -operator  model  _ :  .5  cm,  - :  1.0  cm,  ****; 

2.0  cm,  (c)  dependence  on  surface  autocorrelation 
lengths,  standard  model,  (d)  dependence  on  surface 
autocorrelation  lengths,  slope -operator  model.  Cross 

directed  length  fixed  at  20  cm,  _ ;  10  cm,  - '.  20 

cm,  ****:  40  cm  for  the  parallel  autocorrelation 
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27.  Coherence.  Dependence  on  aucocorrelacion  lengths 

_ :  20,20  and  40,20.  - :  20,40  cm  for  cross- 

directed  and  parallel  autocorrelation  lengths, 

(a)  standard,  (b)  slope-operator 


CHAPTER  V 


IN-PLANE  SCATTER  USING  THE  PIERSON-STACY  OCEAN  SPECTRUM  (42,43) 

The  aim  is  to  illustrate  and  compare  the  far- field  results 
obtained  from  the  models  developed  in  Chapter  3  when  the  surface 
involved  is  the  fully  developed,  non- directional  ocean  spectrum  -  see 
Appendix.  For  convience  we  list  the  most  salient  changes  in  the 
variables  owing  to  the  geometry  (see  Figure  2)  and  the  symmetry  in 
the  spectra; 


l+sintfsin^i 
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l+sintfsinS2 
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(5.3) 


and  the  spectrum  W(x)  is  equal  to  S(/c)//c  for  the  non-directional 


case . 


A.  GENERAL  CASE 
1 .  Standard  Model 

a.  Sllfhtlv  Rough  Ocean 

Implementing  the  above  changes  in  equations  (3.25)  and  (3.26), 
the  covariance  and  intensity  for  in-plane  scatter  are 
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To  evaluate  the  integrals  with  an  omnidirectional  spectrum,  we 


perform  the  change  of  variables 


-  KCOSCl 


Xx  ”  Ksinfl 


and  proceed  using  polar  coordinates.  The  integral  is  now 


-  ^  k2/92v2  CO  2ir  -  ^k2u2  -  ^[k2(V2-U2)cos2n+2k/9/cV2cosn] 

e  ^  J/cdk  f  dnU(K,n)e  e 


(5.6) 


Where  the  omnidirectional  spectrum  S(k)/k  and  the  angular  integration 
is  performed  using  a  procedure  outlined  by  Beckmann  (43,  p.  195). 

The  expression  transforms  into 


-  ^k2^2V2  =0  (r+  ^) 
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When  ^  is  zero,  the  expression  takes  the  form 
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The  evaluation  of  (5.4)  and  (5.5)  with  (5.7)  needs  to  be  performed 
numerically. 


b .  Very  Rough  Ocean 


For  covariance  eq.  (3.30)  is 


A2k2fi(0)f2(e)e 


ik(ri-r2)  -  ^[V2(^2+4^2)/b  )]  . 
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and  the  intensity  is  thus 


A^k^f^re^e  ^ 

4ir2(r'r)2BiB2 


kiM' 

2  UaJ 


(5.9) 


Bi  -  yi-GTiU2 


Bj  -  yi-GTjV^ 


the  square  roots  are  to  be  taken  positive.  For  an  omnidirectional 


spectrum  Ti-T2-T/y2  ,  where 
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A  source  is  placed  500  cm  from  the  ensonified  surface  and  is 
operating  with  a  3°  beamwidth.  The  windspeed  measured  at  1000  cm  is 
514.4  cm/s  (10  knots).  This  10  knot  wind  produces  a  surface  with 
rms  height  of  15.8  cm  and  -7.1*10*®  for  the  parameter  T,  eq.  (5.10). 
Under  these  conditions  the  scattering  strength  varies  with 
source/receiver  polar  angle  for  the  specular  test  as  shown  in  Figure 
28a.  Under  the  same  conditions  the  backscatter  case  is  illustrated 
in  Figure  28b  and  for  the  forward  scatter  case,  28c.  Comparison  of 
the  two  plots  gives  an  idea  of  the  greater  significance  of  the 
coherent  contribution  as  the  frequency  is  increased.  They  also 
display  the  phenomenon  of  a  decrease  in  amplitude  and  broadening  of 
the  curves  as  the  frequency  is  decreased.  Note  as  well  that  the 
^lmplitude  shifts  to  the  specular  direction  -  45°  -  as  the  frequency 
is  increased  in  the  forward  scatter  test. 

If,  on  the  other  hand,  when  the  frequency  and  beamwidth  are  10 
KHz  and  3° ,  respectively,  and  the  measured  windspeed  at  1000  cm  is 
257.2  cm/s ,  514.4  cm/s,  771.6  cm/s  (corresponding  to  5,  10,  and  15 
knots)  the  scattering  strength  will  vary  as  illustrated  in  Figures 
29a,  b,  c.  The  primary  effects  are  a  broadening  and  decreasing 
maximum  amplitude  with  increasing  rms  height  and  slope.  The  forward 
and  backscatter  curves  show  the  fundamental  shortcoming  of  the 


Standard  Model,  namely,  the  severe  misprediction  for  low  grazing 
angles . 


2 .  Coherence 

Figure  30  shows  the  variation  in  the  vertical  coherence  as  a 
function  of  windspeed.  the  source  and  a  fixed  probe  are  200  cm  away 
at  a  depression  angle  of  45“ .  The  receiver  changes  in  depth  from  50 
cm  to  300  cm.  As  would  be  expected,  the  coherence  narrows  with  an 
increase  in  the  surface  rms  slope  and  height.  The  coherence  tends  to 
narrow  as  well  when  the  frequency  or  beamwidth  is  increased  (see 
Figures  31a,  b) .  The  nature  of  this  phenomenon  has  been  discussed  in 
the  previous  chapter. 

B.  SPECIAL  CASE;  WHEN  THE  SPECTRAL  INTEGRALS  CONTAIN  A  GENERALIZED 
FUNCTION 
The  integral 

J  dx^  /  dv  Ne  ^  y  -1  (5.11) 

A  y 

-00  -  00  '' 

satifies  the  equality  when  and  are  very  large  and  when  N  is 
equal  to  A/2ff^ ,  because  in  the  limit  of  very  large  U*  and  the 
integrand  is  a  delta  function.  This  fact  can  be  used  to  explicitly 
evaluate  some  cases  that  involve  an  arbitrary  spectrum. 

In  this  section  the  above  fact  is  used  to  solve  the  equations, 
developed  in  the  previous  chapter,  which  involve  a  non-gaussian 
spectrum.  The  Pierson-Stacy  spectrum  is  chosen  here  because  it  is  an 
interesting  case  involving  a  non-gaussian  spectrum  belonging  to  an 


Vertical  coherence  in  the  backscatter  direction. 

variation  with  windspeed.  _ :  2.57  cm/s ,  - :  514.4 

cm/s,  ****:  721.6  cm/s  measured  at  1000  cm.  Source  at 
200  cm  at  45°  with  3°  beamwidth  and  operating 
frequency  of  10  KHz. 


122 


assumed  Gaussian  height-distributed  Dirichlet  surface.  The  in-plane 
case  is  considered,  as  opposed  to  the  general  geometry,  because  the 
amount  of  algebraic  manipulation  is  reduced  considerably.  However, 
the  approximations  to  be  made  in  this  section  do  not  rely  on  the  fact 
that  the  source -surface  and  surface -receiver  rays  be  co- planar. 

The  integrand  in  eq.  (5.11)  is  not  present  in  all  three  models 
considered  in  the  last  chapter  in  the  form  presented  above.  Thus, 
the  constant  N  shall  be  different  for  all  three  models,  depending  on 
the  other  parameters  of  the  quadratic  term  in  the  exponential  that 
behaves  delta- like  in  the  above  approximation.  The  parameters  and 


are  contained  in  all  three  models  and  one  must  be  careful  to 
consider  cases  of  them  large  not  by  making  the  aperture  wide  but 
rather,  by  making  the  values  of  the  source -surface  distance  very 
large  in  comparison  to  U  and  V  (see  eq.  2.7)  while  maintaining  a  very 
narrow  beamwidth.  This  is  because  an  assumption  common  to  all  three 
models  is  that  the  ensonified  area  dimensions  must  be  comparatively 
small  -  see  Chapter  2 . 

1.  The  Standard  and  Slooe-Ooerator  Models 


htlv  Rough  Surface 


Using  (3.25)  and  (3.26)  as  well  as  the  integration  properties  of 
generalized  functions,  the  Standard  Model  yields 
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for  the  covariance  and  intensity  respectively.  W  is  the  cartesian- 
surface  spectrum. 

The  Slope-Operator  Model  in  the  Fraunhofer  approximation,  Eq. 
(3.59),  gives 
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and  (3.60)  yields 
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for  covariance  and  intensity.  In  both  equations  the  highest  order 
correction  is  neglected. 

b .  Very  Rough  Surface 

For  the  Standard  Model  the  covariance  and  intensity  are 
respectively 
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using  (3.30)  and  (3.31). 


The  Slope-Operator  Model  expressions,  using  (3.61)  and  (3.62) 


are 
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With  GTj^  and  GT2  very  small.  The  added  constants  come  from  the 
normalization  factor  required  for  this  particular  case. 

2 .  Compos ite-Rouehness  Model 

The  two  alternatives;  using  eq.  (3.68),  which  is  the 
conventional  approach,  leads  to 
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for  covariance  and  for  intensity 
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<Tl  is  the  rms  height  and  a'  and  a'  are  the  nns  slopes  of  the  large - 

X  y 

scale  surface,  u  is  the  Heaviside  function.  The  quantity  RjRz  can 
be  replaced  by  Ga^a^a'^  ,  a  quantity  that  must  remain  small  for  these 
expressions  to  be  valid. 

When  the  Slope-Operator  is  incorporated,  the  integrations  in 
(3.72)  and  (3.73)  yield 
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3 .  Discussion 

The  spectral  part  of  the  above  equations  can  be  recast  in 
cylindrical  polar  coordinates,  which  is  more  convenient  than  cartesian 
when  the  surface  is  omni-directional.  In  that  case  the  following 
changes  occur: 
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dYW(y)  -  xd/cdfl 
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where  rj  —  y/c^+2/ck/9cosfl+k^/8^  ,  the  positive  root. 

The  spectral  integrals  will  span  over  k  from  0  to  infinity  and  fl 
from  0  to  2ir ,  except  in  the  Composite-Roughness  Model,  in  which  the 
span  of  K  is  from  to  infinity. 

Figures  32  a  and  b  are  plots  of  the  backscattered  and  forward 
scattered  scattering  strength  as  per  equation  (5.17).  In  these 
computer  runs  the  source  and  receiver  were  set  at  1000  cm  from  the 
ensonified  surface.  The  wind  1000  cm  above  the  surface  was  514.4  cm/s 
and  the  source  was  given  a  1°  beamwidth.  In  the  forward  scatter  test 
the  source  had  a  45°  depression  angle.  It  is  apparent  in  Figure  b  that 
far  away  from  the  specular  direction,  i.e.,  for  receiver  angles  greater 
than  120°  or  so,  the  model  begins  to  fail  and  the  effect  is  more  severe 
for  lower  frequencies.  The  reason  for  this  is  that  for  small  U  and  V, 
the  source  wavenumber  k  must  be  large  in  order  for  the  conditions  kV»l 
and  kU»l  to  prevail. 

A  comparison  of  the  prediction  of  equation  (5.17)  with  that  of 
(5.19)  in  the  backscatter  test  case  is  shown  in  Figures  33a  and  b. 

These  are,  respectively,  the  Standard  Model  and  the  Slope -Operator 
Model  for  source/receiver  distance  of  1000  cm,  1  degree  beam  and  20 
KHz.  The  wind  is  514.4  cm/s  at  1000  cm  above  the  surface.  The  Slope - 
Operator  Model  and  the  Composite-Roughness  Model  using  equation  (5.21) 
are  compared  next  in  Figures  34a  and  b.  In  this  case  the  source  is  70 
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Figure  32.  Scattering  strength.  Eq.  (5.17)  (a)  backscatter,  (b) 
forward  scatter.  Receiver/source  at  1000  cm.  1° 

beamwidth.  Wind  at  1000  cm  is  514.4  cm/s.  _ :  15 

KHz,  - :  20  KHz,  ****:  50  KHz. 
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Figure  33.  Backscatter.  (a)  Eq .  (5.17)  standard  model,  (b)  eq. 

(5.19)  slope-operator  model,  20  KHz  1°  beamed  source 
Wind  at  1000  cm  is  514.4  cm/s 


cm  away  and  all  other  parameters  discussed  above  are  held  fixed. 

Figure  34b  shows  the  characteristic  Bragg  resonance  dip  in  the 
backscatter  plots,  which  are  typical  of  the  experimental  findings.  The 
dip  is  more  noticeable  when  the  beamwidth  is  narrow  -  see  Figures  35a 
and  b. 

Figures  35a,  b,  and  34b  have  all  been  plotted  using  a  quarter 
source  wavenumber  value  for  spectral  cutoff  value  .  The  choice  of 
cutoff  is  not  at  all  given  by  the  theory  itself.  Figures  36a,  b,  and  c 
show  the  variation  in  the  backscatter  with  choice  of  cutoff  wavenumber. 
In  these  figures  the  frequency  is  10  KHz  and  the  choice  of  is  (a) 
equal  to  the  source  wavenumber,  (b)  .8  times  the  wavenumber,  (c)  0.25 
times  the  wavenumber. 

Another  feature  of  the  Composite -Roughness  Model  as  per  equations 
(5.21)  and  (5.23)  is  that  care  must  be  exercised  in  order  to  not 
operate  out  of  the  bounds  of  the  regime  of  validity.  Assuming  the 
beamwidth  is  small,  the  four  conditions  that  must  be  satisfied  are: 

(a)  GCg  «  1,  which  leads  to  «  1.  This  permits  the 

series  expansion  of  the  exponential  involving  GCg  and  the  retention  of 
the  first  two  terms. 

(b)  cr"^R^  «  24<7^  -  12a' ^R^,  where  R  is  a  radial  distance  on 

L  L  I, 

the  insonified  area.  In  order  for  the  Taylor  series  expansion  of  the 
large-scale  surface  autocorrelation  to  remain  valid,  eq.  (3.66),  the 
terms  higher  than  second  order  must  be  insignificant  compared  to  those 


retained. 
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(c)  U  and  V  must  be  of  comparable  or  smaller  size  than  R.  This 

is  necessary  because  an  early  assumption  of  the  TPM  was  that  the 

integral  over  the  entire  surface  was  the  same  as  the  integral  over  the 

insonified  portion  due  to  the  minute  contribution  of  the  integrand 

outside  of  this  region.  As  well,  it  is  necessary  for  the  quadratic 

terms  in  the  spatial  integrations  of  (3.70)  to  be  such  that  U"^  » 

and  V  »  Ga''^/a^  in  order  to  permit  the  integration  over  all 
X  y 

space . 


(d)  Related  to  the  generalized  function  integration,  it  is 

necessary  that  V  ^ 8 a «  1/2  and  U  ^  +4k2sin^0a^2  «  1^/2 

because  the  delta- like  integration  becomes  possible  when  2R^/V^  «  1 

and  2R^/U2  «1. 

v 

The  Pierson-Stacy  report  on  ocean  spectra  (42)  is  to  the  author's 
knowledge  the  best  study  available  in  the  open  literature.  The  report 
is  not  considered  to  be  a  definitive  study  on  the  subject  owing  to  the 
fact  that  it  is  a  compilation  of  reports  made  by  several  investigators 
who  considered  limited  surface  wavenumber  ranges  in  different 
experimental  conditions,  using  a  wide  range  of  measuring  techniques. 

In  fact,  most  investigators  agree  that  a  thorough  and  reliable  study  on 
the  ocean  spectrum  would  be  welcomed,  but  the  task  has  not  been 
performed  due  to  the  challenges  involved  with  the  measurement  of 
spectra  in  the  open  sea. 

This  chapter  has  been  devoted  mainly  to  the  in-plane  case,  hence 
the  use  of  the  omnidirectional  spectrum  W(x)  -  S(k)/k.  The  non- 
coplanar  case  would  be  most  interesting  and  general,  if  the  directional 

I 

f 
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spectrum  is  utilized  instead  of  the  omnidirectional  one  used  here. 

This  would  require  making  none  of  the  geometrical  simplifications  used 
in  the  in-plane  situation,  as  well  as  the  incorporation  of  angular 
variation  in  the  spectrum  function. 

Angular  variation  in  the  spectrum  would  be  expressed  in  the  form 

V(k.Q)  -  S(x)F(«.fl)  (5.29) 

with  S(x)  being  the  non-directional  function  used  in  this  study  and 
F(/c,fl)  most  researchers  aegree  has  a  general  form 

F(K,n)  a  G(s)|cos^  (n-0')|2»  (5.30) 

where  s  is  a  function  of  k  such  that  F(x,n)  tends  to  a  narrower  and 
higher  amplitude  angular  distribution  (44)  and  G(s)  is  a  normalization 
function  such  that  the  integral  over  all  wavenumber  space  of  equation 
(5.29)  is  one.  The  constant  0'  is  aligned  with  the  mean  wind  velocity 
direction. 

The  surface  autocorrelation  can  be  a  function  of  time  as  well. 
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This  implies  that  the  surface  has  a  well  defined  temporal  dependence 
and  that  the  covariance  will  be  space  as  well  as  time  dependent.  The 
autocorrelation  function  can  have  an  envelope  which  travels  in  time 


with  same  group  velocity.  .A  case  in  point  (38,  p.  518)  would  be  the 


wind-blown  water  waves,  for  which  the  autocorrelation  can  be 


C  -  exp 


cos  n(^  -  r)Jexp(-bx^) 


(5.31) 


A 

where  the  waves  are  moving  in  the  y  direction  with  phase  velocity  v  and 
angular  frequency  0.  The  envelope  travels  with  group  velocity  u  and  is 

A 

coherent  within  an  approximate  range  r  in  time.  In  the  x  direction  it 
dies  off  by  virtue  of  the  parameter  b,  the  length  of  each  crest. 

In  this  type  of  highly  organized  surface  states,  the  spatial  and 
temporal  components  of  the  autocorrelation  are  not  separable.  In  a 
confused  sea,  the  waves  travel  in  all  directions  and  we  would  expect 
the  spatial  correlation  function  to  be  isotropic,  as  was  assumed  in 
this  study,  and  that  the  spatial  and  temporal  components  be  separable. 


CHAPTER  VI 


SUMMARY  AND  SUGGESTIONS  FOR  FURTHER  STUDY 


In  this  study  two  different  philosophies  for  representing  the 
scattered  field  when  a  randomly -rough  surface  is  involved  have  been 
presented.  The  MSP  is  simple,  adequate  when  the  roughness  is  very 
slight.  It  doesn't  rely  on  the  Kirchhoff  approximation  and  will 
tolerate  steep  slopes.  The  scatter  is  diffuse,  whereas  the  TPM 
consider  it  more  of  a  distortion.  The  TPM  does  not  handle  steep  slopes 
but  can  work  well  with  rough  surfaces.  The  drawbacks  of  the  latter  are 
the  use  of  the  Kirchhoff  approximation  as  well  as  cumbersome 
mathematics . 

Within  the  TPM,  several  variants  have  been  presented.  The 
Standard  Model,  simplest  of  all,  is  ill  behaved  for  angles  far  from 
the  specular  direction  because  the  incoherent  component  makes  bad 
estimates  of  its  contribution  and  because  an  assumption  in  the 
integration-by- parts  techniques  is  that  the  angular  dependence  is  a 
function  which  is  nearly  constant  so  that  that  it  can  be  removed  from 
the  integral. 

The  other  extreme  in  complexity  is  the  Slope-Operator  Model,  which 
is  shown  to  behave  properly  for  all  grazing  angles.  This  model  is 
compared  when  the  surface  is  slightly  rough  to  the  case  when  the 


gradient  of  the  normal  of  the  surface  is  approximated  by  d/dz.  In 
effect  this  is  a  comparison  between  the  full  model  and  the  coherent 
component  itself  and  the  two  obvious  shortcomings  are  that  the  latter 
does  not  incorporate  any  spectral  information,  nor  does  it  apply  to 
other  than  low  frequency  cases.  The  Slope-Operator  model  is  presented 
using  both  the  Fraunhofer  and  Fresnel  phase  approximation.  The  Fresnel 
case  is  more  complicated  and  presents  no  more  information  than  the 
Fraunhofer  model  in  the  farfield.  In  principle  the  Fraunhofer  case 
will  diverge  in  the  specular  direction  if  no  source  beampattern  is 
used.  This  is  because  the  integral  over  the  ensonified  area  is  assumed 
to  be  the  same  as  the  integral  between  minus  and  plus  infinity,  since 
the  integrand  with  a  Gaussian  beampattern  restricts  significant 
contributions  to  the  ensonified  area.  The  Fresnel  approximation 
presented  in  this  study  is  derived  using  the  binomial  expansion.  This 
alternative  is  admittedly  poor  when  compared  to  a  power  series 
expansion  in  the  phase  but  is  considerably  easier  to  implement. 

Lastly,  the  MSP  and  TPM  are  effectively  merged  in  the  Composite- 
Roughness  Model.  It  is  rather  easy  to  implement  and  offers  good 
results  when  the  surface  spectra  can  be  divided  into  two  distinct 
regimes.  However,  the  model  is  very  sensitive  to  the  choice  of 
spectral  partitioning  when  it  is  not  apparent  where  the  division  in  the 
spectrum  is  to  be  made. 

Expressions  for  the  covariance  and  intensity  were  developed.  The 
coherent  component  of  the  intensity  is  shown  to  be  a  good  estimator  of 


surface  roughness.  Spectral  content  turns  up  in  the  incoherent 
component.  The  covariance  of  the  field  is  shown  to  be  analogous  to  the 
mutual  covariance  quantity  whereas  the  intensity  is  the  same  as  the 
self -covariance . 

No  attempt  was  made  to  conduct  an  in-depth  coherence  study,  but 
several  examples  were  provided  so  that  the  models  can  be  compared  and 
some  qualitative  features  of  coherence  calculations  can  be  highlighted. 
The  coherence  is  shown  to  drop  with  increasing  beamwidth  and  frequency. 
This  feature  is  useful  in  interferometric  partial  coherence  studies. 
When  the  Fraunhofer  phase  approximation  is  used,  the  coherence  phase  is 
determined  solely  by  the  wavenumber  times  the  difference  in  magnitude 
of  the  distances  between  the  surface  and  the  two  receivers.  The  point 
was  made  that  the  horizontal  coherence  tends  to  vary  very  little 
compared  to  the  vertical  coherence  because  the  random  processes 
considered  are  a  function  of  vertical  coordinates.  Coherence  is  a  good 
estimator  of  source  localization  and  it  provides  considerable  spectral 
information. 

Although  partial  coherence  is  a  well-developed  theory  in  the  field 
of  optics,  in  acoustics,  a  comprehensive  theory  encompassing  the 
fundamental  notions  is  yet  to  be  formulated.  The  acoustics  community 
would  undoubtedly  welcome  a  monograph  in  this  area. 
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In  optical  coherence  studies,  rapidity  of  phase  fluctuations  is 
the  primary  challenge,  whereas  in  acoustics  the  difficulty  is  that  near 
perfect  shadows  are  unattainable  in  an  important  frequency  range. 

These  aspects,  along  with  the  distinctly  different  near-field 
characteristics  and  the  polarization,  must  be  dealt  with  differently  in 
acoustics . 

Reliable  experimental  data  on  coherence  measurements  involving 
Gaussian  surfaces  as  well  as  the  ocean  surface  is  needed  in  order  to 
make  an  assessment  of  the  validity  of  the  several  models  presented  in 
this  study.  Exact  numerical  evaluation  of  equations  (2.1)  and  (2.2) 
would  permit  an  analysis  of  the  validity  and  interaction  of  the  many 
assumptions  and  approximations  made  in  this  study,  as  well  as  a  means 
of  comparing  the  two  equations  for  similarity  in  results.  A  difficult 
but  necessary  project  is  the  measurement  of  the  ocean  spectra  under 
more  controlled  and  consistent  conditions.  The  tabulation  of  the 
surface  autocorrelation  could  be  performed  as  once  did  Fortuin  and 
deBoer  (12) ,  Latta  (45) ,  and  Mellen  (46)  with  the  outdated  Roll-Fisher 
and  Neumann- Pierson  spectra. 

Time  dependence  can  be  introduced  in  the  problem,  so  that  a  model 
for  spatial  and  temporal  covariance  can  be  formulated.  Time  coherence 
would  enable  the  determination  of  finite  spectral  width  of  the 
scattering  surface  and  surface  movement,  in  addition  to  surface 
autocorrelation  "movement."  The  directional  spectrum  and  full 
geometric  implementation  would  be  most  useful  in  this  case. 
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A  logical  and  extremely  useful  continuance  of  this  work  would 
be;  (1)  the  quantitative  analysis  of  the  sensitivity  of  the  models 
to  field 

fluctuations,  (2)  the  importance  of  higher-order  multiple  scatter 
terms  to  coherence  calculations,  (3)  the  contribution  of  multiple 
scatter  to  field  flucations,  and  (4)  the  development  of  covariant 
expressions  that  are  valid  in  the  near- field. 

Derivation  of  expressions  for  higher-order  moments  (which  are 
less  sensitive  to  inhomogeneities  in  the  media)  is  required  for  a 
theory  of  partial  coherence,  particularly  if  the  general  theory 
for  Gaussian  as  well  as  non-Gaussian  surfaces  is  formulated. 

Though  higher  order  moments  do  not  provide  any  new  information  if 
the  surface  is  Gaussian,  they  do  so,  however,  for  non-Gaussian 
surfaces.  The  theory  underlying  non-Gaussian  surface  statistics 
can  be  found  in  Beckmann's  paper  (38). 
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APPENDIX 


THE  PIERSON- STACY  SPECTRUM 


The  non- directional  ocean  wavenumber  elevation  spectrum  is 
included  here  for  convenience.  This  summary  has  been  extracted  word 
for  word  from  a  summary  prepared  by  McCammon  (48)  on  the  subject  of  the 
Pierson-Stacy  (42)  Reports. 

This  wavenximber  elevation  spectra  is  non-directional  with  the 
units  of  cm^  and  is  related  to  the  RMS  waveheight  by 


«T^>  -  /  W(K)dK. 


The  spectral  equations  for  each  region  are 
(1)  Pierson  Moskowitz  region 

Wi(K)  -  o  exp[-B/K2)/(2K3). 


(2)  Stacy  region 


(A.l) 


(A. 2) 


Ws(K)  -  (2.04xl0®+1.02xl0*U*)exp[-79.2K^'^]/K^^^  (A. 3) 


(3)  Kitaigorodskii  region 


W2(K)  -  a/(2^''^K’'^). 


(4)  Leykin- Rosenberg  range 


W3(K)  -  W,  (K3)(K/K3)‘^. 


(A. 4) 


(A. 5) 
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.  1  Pierson  wavenumber  elevation  spectr'-im 
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(5)  Mi tsuyasu- Honda  range 


W*(K)  -  0.875(2ir)'’'^(l+3TK2)/(g''’‘^’^^(K+K3r)'‘’‘"^’'^  ). 

(A. 6) 

(6)  Cox  viscous  cutoff  range 

WjCK)  -  7.04x10'^  U**/K».  (A. 7) 
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where  a 


q 

p 


8.1x10 

0.74  g2/U<  , 

51.71A^*^ 

5.0  -  log^^[U*). 


f  -  7.58x10 

g  -  980. 


These  regions  are  linked  together  in  K  space  as 


W(K)- 


max( 


max( 


mini 


Wi(K) 

Ws(K) 

Wi(K) 

Ws(K) 

W2(K) 

WjCK) 

W*(K) 

WjCK) 


U*  >  35.8  else  use 
Wi(K) 


U*  >  75.76  else  use 
W^CK) 


0<K<Ki 

Ki<K<Ki 

Ki<K<K2 

K2  <K<K3 
K3>K 


.  » 


-  g- 


'>./i 


w*.wv: 
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Ki  -  0.702g/U^g^ 

Ki  -  51.71/U*2, 

Kj  -  0.359, 

K3  -  0.942, 

U  -  wind  speed  at  19.5  m, 

U*  -  friction  velocity  obtained  from  the  equations 

U2  -  U*  ln(z.Zo)/0.4, 


Zq  -  0.684/U*  +  4.28xlO‘^U*2  .  0.0443. 


The  units  of  these  equations  are  CGS,  W(K)  is  given  in  cm^ , 


cm  ^ ,  U  in  cm/sec , 


pjL  m 

l>r/c. 
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